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Abstract. We apply ideas from C*-algebra to the study of disordered topological insu- 
lators. We extract certain almost commuting matrices from the free Fermi Hamiltonian, 
describing band projected coordinate matrices. By considering topological obstructions to 
approximating these matrices by exactly commuting matrices, we are able to compute in- 
variants quantifying different topological phases. We generalize previous two dimensional 
results to higher dimensions; we give a general expression for the topological invariants for 
arbitrary dimension and several symmetry classes, including chiral symmetry classes, and 
we present a detailed -ftf-theory treatment of this expression for time reversal invariant three 
dimensional systems. We can use these results to show non-existence of localized Wannier 
functions for these systems. 

We use this approach to calculate the index for time-reversal invariant systems with spin- 
orbit scattering in three dimensions, on sizes up to 12 3 , averaging over a large number of 
samples. The results show an interesting separation between the localization transition and 
the point at which the average index (which can be viewed as an "order parameter" for the 
topological insulator) begins to fluctuate from sample too sample, implying the existence of 
an unsuspected quantum phase transition separating two different delocalized phases in this 
system. One of the particular advantages of the C*-algebraic technique that we present is 
that it is significantly faster in practice than other methods of computing the index, allowing 
the study of larger systems. In this paper, we present a detailed discussion of numerical 
implementation of our method. 
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1. Topological Insulators and Wannier Functions 





Consider a free Fermi Hamiltonian, described by a matrix "H^-, where i,j label sites on 
some lattice. Suppose the Hamiltonian H is local, so that 7Uj decays rapidly in the spacing 
between i and j, and that H has a gap in its spectrum. Then, the system can be in 
either topologically trivial or topologically nontrivial phases. An example of a topologically 
nontrivial phase is provided by the quantum Hall effect in two dimensions; if the Hamiltonian 
if is a lattice Hamiltonian in a magnetic field with a non-zero Hall conductance, then there is 
a topological obstruction to continuing the Hamiltonian if to a trivial Hamiltonian without 
either making the gap small or violating locality. Here, by a "trivial Hamiltonian", we mean a 
Hamiltonian in which all sites are decoupled, so that T-tij is diagonal, while the magnitude of 
the gap is related to the system size later. Equivalently, such an obstruction is an obstruction 
to having localized Wannier functions in the system[T], as discussed further below. 

We can consider also such a system with symmetries imposed. For example, we can require 
that the Hamiltonian have time reversal invariance. Such time reversal invariant topological 
insulators were considered in [5]. Such a system is topologically trivial when considered as a 
time reversal non-invariant system; that is, there is no obstruction to continuing the system 
to a trivial system if we do not require that the path connecting the Hamiltonian to a trivial 
Hamiltonian also be time reversal invariant. However, there is an obstruction to continuing 
such systems to time reversal invariant trivial systems along time reversal invariant paths. 

Many other symmetry classes can be considered, and in addition one can consider systems 
in higher dimensions [6]. One important concept in the classification of topological systems in 
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higher dimensions is the idea of "stable equivalence". We consider two Hamiltonians "Ho>%i 
to be connected if we can add some number of additional trivial degrees of freedom to H®, 
then follow a continuous path in parameter space maintaining locality and spectral gap, 
and then arrive at a final Hamiltonian which is equal to "Hi plus possibly additional trivial 
degrees of freedom. Using this definition of stable equivalence, we define a system to be 
in a topologically nontrivial phase if it is not stably equivalent to a trivial Hamiltonian. A 
general table of topological obstructions in different symmetry classes and dimensions was 
presented in [8] . 

In [TUHJUD], an alternative approach to studying these topological insulators was presented, 
based on C*-algebras. We now describe the approach for systems in the three classical uni- 
versality classes, where the Hamiltonian is either an arbitrary Hermitian matrix, a real sym- 
metric matrix, or a self-dual matrix, respectively. These classes are sometimes referred to as 
A, AI, AIL We use the names GUE,GOE,GSE (Gaussian unitary, orthogonal, and symplec- 
tic) for these classes; this is not intended to imply that the Hamiltonian is drawn from some 
particular Gaussian distribution but is simply shorthand for the time-reversal symmetries 
imposed on the Hamiltonian. Later in this paper we give the appropriate generalization for 
the chiral cases (which we refer to as chiral, chiral real, and chiral-self-dual) but we begin 
with the three classical classes. First, consider a system on a d- dimensional sphere (the case 
of a system on a torus or on other topologies is presented in a later section of the paper; the 
sphere case is the simplest to explain first). That is, we imagine each lattice site i as located 
somewhere on the surface of a d- dimensional sphere, with coordinate xi(i),x 2 (i), ...,Xa+i(i)- 
We normalize the radius of the sphere to unity, so that ^2 a x a (i) 2 = 1. Let Xi, X d+ i be 
coordinate matrices. These are diagonal matrices, with diagonal entries {X a )n = x a (i). In 
the GSE case, there are two different spin-states per site, and so these coordinate matrices 
X a are self-dual matrices. Thus, 



where / is the identity matrix. We then compute the projector P onto the space of occupied 
states. This is the space of eigenvalues of the Hamiltonian with energy less than the Fermi 
energy, Ep. We will see later that the properties of the system will depend on the value of Ep 
chosen, and we will observe phase transitions as a function of Ep, as the system changes from 
an ordinary insulator to a diffusive metal, to a topological insulator (experimental variation 
of the Fermi energy is possible via gating in some such systems [3T]). 

We then define a set of band projected matrices, H r . Let PX a P be the "band projected 
position matrix". We write 



where the two blocks in the above matrix correspond to the space spanned by the kernel of P 
and the range of P (the space of empty states and occupied states, respectively). Thus, the 
size of the matrix H r is equal to the number of occupied states, which will be important in 
reducing the numerical effort later. Note that if H is GUE,GOE, or GSE, then the operators 
H r are also GUE, GOE, or GSE respectively, so that they inherit the symmetries of H. 

Suppose the commutators [P, X a ] are small. This occurs if the Fermi energy E F is in a 
spectral gap or in a mobility gap. In particular, to determine the size of the gap needed 



d+l 





a=l 
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to make [P, X a ] sufficiently small, we have to consider the ratio between the range of the 
Hamiltonian and the system size. Above we have normalized distances so that 'Yl ia x a {i) 2 — 1 
and the sphere has radius unity. In this case, if there are many sites on the sphere, then the 
distance between each site is small, tending to zero as the number of sites tends to infinity. 
For clarity in the present discussion, we prefer a more general normalization of distances, so 
that we can instead normalize the distance between sites to a constant, independent of the 
number of sites, and change the linear size of the system together with the number of sites. 
To do this, let us instead normalize so that 'Y^ a x a {i) 2 = L 2 '', for some length scale L, and 
let H r = PX a P/L, so that the previous discussion corresponds to the choice L — 1. Let 
us assume we have finite range interactions, so that Hij vanishes if the distance between i 
and j is larger than some interaction range R, which is held fixed for all system sizes (our 
numerical studies below correspond to a choice R = 1 since they involved nearest neighbor 
hopping, with the distance between sites normalized to 1). Let us assume that \\H\\ < J for 
some constant J. Then, if Ep is is separated by a gap AE from the spectrum of H, one can 
show that the commutator \\[P, X a ]\\ is bounded by a constant times RJ/(LAE). 

Since [P, X a ] is bounded, the band projected position matrices almost commute with each 
other and almost square to the identity: 

(1.3) \\[H r ,H s }\\^0, 

(1.4) 5» J - 

r 

We refer to this as a "soft" representation of the sphere S d . We sometimes quantify the 
approximation in the above equations, saying that a set of matrices form a ^-representation 
of the sphere if 

(1.5) \\[Hr,H s ]\\ <8, 

(1.6) llE^ 2 - J ll^- 

r 

One can show that we have a 5-representation of the sphere with 5 of order jl] 

(1.7) 5 < const, x (RJ/LAE) 2 

= const, x (vlr/LAE) 2 , 

where the quantity vlr has units of velocity (we chose to introduce this velocity because 
this allows one to also describe Hamiltonians which do not have finite range interactions but 
instead have interactions that decay sufficiently rapidly with distance; for such Hamiltonians 
one can prove the same bound with only a small amount of extra work). 

A well-studied problem in C*-algebra is whether such a set of almost commuting matrices 
can be approximated by a set of exactly commuting matrices. In the absence of symmetry, 
this problem is completely understood in the case of the two-sphere and the two-torus. The 
answer is that the approximation is possible if and only if a certain topological invariant, 
discussed in the next section, vanishes. This topological invariant is an integer, and may be 
(see the section on mathematical problems for more detail on this question) identified with 
the Hall conductance. In other symmetry classes and other dimensions, we have identified 
other topological invariants, presented in the next section. In one particular case (the case 
of fermionic systems with charge conservation but no other symmetries, which corresponds 
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to the case in which the X a are arbitrary Hermitian matrices with no other symmetries), it 
is possible to determine the complete set of topological invariants of matrices in the stable 



Our main claim is that the topological invariants of these almost commuting matrices can 
in general be identified with the topological invariants of free Fermi system; we have not 
proven this in all cases, but we have observed several cases fitting in the general pattern and 
have a proof in some cases of this identification. 

In this paper we begin by describing the relationship between the topological invariants of 
the almost commuting matrices and the existence of Wannier functions in the free fermion 
problem, and sketch the outline of our numerical procedure. We then describe how to 
compute invariants of the matrices for a spherical geometry in several cases using the so- 
called "Bott matrix". We summarize various mathematical questions, and then discuss how 
to compute invariants in other geometries. We then describe chiral classes, where invariants 
of almost commuting matrices can again be used to classify topological phases, but where the 
procedure of constructing the almost commuting matrices is different, so that the relevant 
matrices are not the band projected position matrices. We then present numerical results on 
a time-reversal invariant topological insulator with disorder in three dimensions. Finally, we 
give additional mathematical details on the K-theory to compute invariants in more general 
cases, and we describe our numerical implementation in detail. 

1.1. Wannier Functions and Trivial Hamiltonians. The topological classification of 
free Fermi Hamiltonians we are interested in is a stable classification as discussed above. The 
question to ask is whether a given gapped Hamiltonian Ho can be connected by a continuous 
path of gapped local Hamiltonians to the trivial Hamiltonian, where the trivial Hamiltonian 
is a diagonal Hamiltonian so that in the trivial Hamiltonian each site is either occupied or 
empty. In the GSE case, we should instead have two states per site, corresponding to spin up 
and down, with the two states having the same energy. Then, such a diagonal Hamiltonian 
is a member of the appropriate symmetry class, either GUE,GOE, or GSE. 

We now show that this is equivalent to the question of whether "Ho has localized Wannier 
functions, where by localized Wannier functions we mean that we can find an orthonormal 
set of functions which are localized in space and which span the range of the projector P 
onto the occupied states of "H - 111 the case of the GOE or GSE we require that the Wannier 
functions respect the symmetry of the Hamiltonian. Thus, the Wannier functions must be 
real in the GOE case and the Wannier functions must occur in time-reversal invariant pairs 
in the GSE case. 

Suppose a given Hamiltonian can be connected to a trivial Hamiltonian %i by a smooth 
path of Hamiltonians H s . Clearly, the trivial Hamiltonian has localized Wannier functions: 
for each site on which the diagonal entry of H.\ is negative, corresponding to an occupied 
site, we have one Wannier function localized on that site. Let this set of Wannier functions 
be {f a (l)}, where a is a discrete index labeling the different Wannier functions. We now 
show that it is possible to construct a set of Wannier functions {v a {s)} for all Hamiltonians 
along the path, using quasi-adiabatic continuation to continue the set of Wannier functions 
along the path [10]. Define 



where T> s is a quasi-adiabatic continuation operator. This gives a set of Wannier functions for 
all Hamiltonians along the path. With an appropriate choice of quasi-adiabatic continuation 



limit [12J. 
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operator, these Wannier functions are superpolynomially localized: for all s G [0,1], the 
amplitude of a given Wannier functions decays superpolynomially away from the site on 
which the function is localized at s = 1. 

Thus, if a Hamiltonian is connected to the trivial Hamiltonian, it has Wannier functions, 
so an obstruction to finding localized Wannier functions implies an obstruction to continuing 
to the trivial Hamiltonian. Conversely, if a Hamiltonian "H has a set of localized Wannier 
functions, {t> a (0)}, we can continue this Hamiltonian to the trivial Hamiltonian as follows. 
First, recalling that we are only interested in stable equivalence, add additional sites to the 
system, one for each Wannier function (so that the number of added sites is equal to the 
number of occupied states). Let these sites be decoupled from each other and the rest of the 
Hamiltonian, with an energy +1 for each site, so that we consider the Hamiltonian 

(1.9) n ®i = 

where the dimension of the second block is equal to the number of added sites. Let i a denote 
the added site corresponding to a given Wannier function v a , and let \i a ) denote the basis 
vector corresponding to this site. Then smoothly continue the first block from Hq to the 
spectrally flattened Hamiltonian 1 — 2P. Then, note that 

(1.10) P = J2K)(Va\- 

a 

Thus, the spectrally Hamiltonian is equal to 1 — 2 J2 a \ v a)(v a \- Follow the continuous path 
of Hamiltonians 

(1.11) l-2j2\Xa(0))(x a (d)\, 

a 

where 

(1.12) \x a (9)) = cos(9)\v a )+sm(9)\i a ). 

from 6 = to 6 = n/2. This gives a continuous path of gapped Hamiltonians to a trivial 
Hamiltonian, and because the Wannier functions are localized, the Hamiltonians are local 
throughout (interaction terms in the Hamiltonian decay superpolynomially if the Wannier 
functions decay superpolynomially). 

Again, this emphasis on Wannier function is specific to these three classical ensembles. 
Later we consider chiral ensembles. In p], it is shown that if localized Wannier functions 
exist, then it is possible to approximate the matrices H r by exactly commuting Hermitian 
matrices. Conversely, given the ability to approximate the H r by exactly commuting Hermit- 
ian matrices, one can define a set of Wannier functions, albeit ones which may be localized 
in a weaker sense than exponential or even than superpolynomial. 

There are topological invariants that can be associated to soft-representations of the 
sphere. By calculating these, we have a mechanism to prove a Hamiltonian H cannot be 
deformed, via a path that stays gapped and local, to a trivial Hamiltonian. The chain of 
implication used in this are the contrapositives of the down arrows in figure 11.11 The invari- 
ants that we compute of the soft sphere will be invariant under any continuous deformation 
of the matrices H r so long as the matrices H r continue to prove a ^-representation of the 
sphere for sufficiently small 5 (5 less than some numeric constant which depends upon di- 
mension). Using the relation between the energy gap and S, this proves that if a Hamiltonian 
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The Hamiltonian H. can deformed via 
gapped and local Hamiltonians to a triv- 
ial Hamiltonian. , 

11(a) 

V II 

(Localized Wannier functions exist. } 

H(b) 

1} II 

The almost commuting Hermitian ma- 
trices Hi, ... , H d+ i are close to exactly 
commuting Hermitian matrices. 

II (c) 

JL ii 

The JT-theory associated to 
Hi, ... , H d+ i equals the K-theory 
of some triple of exactly commuting 
v Hermitian matrices. y 

Figure 1.1. We indicate the directions of implication between various dif- 
ferent mathematical concepts. The bold arrows pointing downward indicate 
directions of implication that we prove to be true as stated. The dashed ar- 
rows indicate directions of implication that are either not proven to be true as 
stated (but a slightly modified form is proven to be true), or are conjecture. In 
particular, a) we prove that the existence of Wannier functions implies stable 
equivalence of Hamiltonians, a weaker statement, in subsection ( 11.11) . (b) In 
PP, we prove that the ability to approximate by exactly commuting matrices 
implies some weaker statements about the localization properties of Wannier 
functions. In particular, we showed that this implies that there exist Wannier 
functions whose variance is small compared to the system size, but we did not 
prove the existence of exponentially localized Wannier functions, (c) This is 
proven (see pQ and references there) in the GUE case for d = 2. It is con- 
jectured in the GOE/GSE case. For the case of more than 3 matrices (that 
is, for systems in more than 2 dimension), we can only hope to have stable 
approximation. See problem ( 13 .4p . 



is continuously deformed to a trivial Hamiltonian then the gap must become of order vlr/L 
somewhere along the path. 

That is, we prove the existence of an obstruction to deforming some Hamiltonian to a 
trivial Hamiltonian by computing some invariant of the matrices H r . This forms the basis of 
our numerical algorithm: we consider a given Hamiltonian and choose a Fermi energy Ep. We 
compute the projector P onto states with energy less than Ep using standard techniques of 
linear algebra (in section (Q we discuss how using symmetries such as time reversal symmetry 
can improve the numerical accuracy of this calculation). We then compute invariants of 
those matrices; these invariants are elements of a group, either Z or Z 2 , describing different 
invariants of the original free fermion problem. 
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The absence of localized Wannier functions for systems in topologically nontrivial phases 
represents a potential difficulty for numerical simulation algorithms based on Wannier func- 
tions such as those in [2] . It has also interesting implications for the application of MERA 
(multi-scale entanglement renormalization ansatz) schemes such as in [3] based on a hier- 
archical construction of wavefunctions for free fermion problems. It is claimed in [3] that 
for gapped systems after a few rounds of MERA for a gapped system, one converges to a 
fixed point describing a product state. However, such a fixed point corresponds to a system 
which does have localized Wannier functions. Thus, when MERA is applied to a topolog- 
ically insulator, one cannot converge to a product state and must instead keep nontrivial 
entanglement at all length scales; converges to such a product state fixed point happens only 
for topologically trivial insulators. We will present this elsewhere, showing that instead one 
gets a structure similar to [4], where a constant number of degrees of freedom must be kept 
at each scale. 

1.2. Real, complex and self-dual matrices. In the GOE and GSE cases, our method 
requires and benefits from preserving the needed symmetry throughout the calculations. It 
is limiting in the GOE case to work only with real matrices so at times we use complex 
symmetric matrices. For example, we will consider a complex matrix that is symmetric and 
unitary, so U T = U and W = U~ l . Mathematical readers should note that U T denotes just 
the transpose, U the conjugate, and W = U the adjoint. When we extract the associate 
sin/cosine pair those matrices will not just be Hermitian, but real symmetric. 

In the GSE case, we need to specify the dual operation on matrices, and make precise the 
meaning of "time reversal invariant pairs" of vectors. The dual operation on M^a^C) is only 
determined up to a choice of a specified matrix Z with the properties 

z- 1 = z f = Z T = -z, 

which means it is real orthogonal with eigenvalues ±i. Unless specified otherwise, what we 
are using is Z = ( ° ~J J . The dual X" of a matrix is then 

(1.13) X s = -ZX T Z. 

We use the operation (j to denote the dual to use a symbol distinct from overline, dagger, 
star, and so on, which we use for other purposes. In terms of N-by-N blocks, 



" A 


B ' 




C 


D 





We need the associated conjugate linear operation T : C 2N — > C 2N defined by 

(1.14) Tv = -Zv. 

For v any vector, Tv is orthogonal to v. 
A matrix W is symplectic if W T ZW = Z. 

Lemma 1.1. Suppose U us a unitary matrix. The following are equivalent: 

(1) U is symplectic; 

(2) W_ = 0; 

(3) ZU = UZ; 

(4) If v is column j of U for j < N then column N + j of U is Tv. 
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Lemma 1.2. If U is a symplectic unitary and X and U are in M 2A r(C) then 

(rfxuY = U ] X^U. 

For low dimensional GSE systems, we derive real matrices whose invariants are equivalent 
to invariants of self-dual matrices. In dimension 4 for a GOE system, we do a conversion 
the other way, so end up with self-dual matrices that encode the invariants of real matrices. 
This is part of how we take advantage of Bott-periodicity in TT-theory, avoiding homotopy 
calculations involving 4-spheres or 5-spheres. 

The trick we use is classic. Tensoring together two dual operations leads to an operation 
that is equivalent to the transpose. This has a simple physical interpretation: the self-dual 
operation is the time reversal operation on a system with half-odd-integer spin. Tensoring 
together two systems with half-odd-integer spin gives a system with integer spin, and the 
time reversal operation on such a system is equivalent to the transpose, after an appropriate 
change of basis. The next technical lemma specifies the appropriate basis change to make 
the time reversal into the transpose. 

Lemma 1.3. Consider 

U = —=(I -iZ <g> Z') . 
v2 

For all X G M 2JV (C7) and Y G M 2m (C), 

C/ t (X a <g> y") U = (C/ f (X <g> Y) u) T 

where (A ® B) T = A T ® B T , i.e. is the usual transpose of AN-by-AN matrices. Here Z and 
Z' are the matrices, of sizes 2N-by-2N and 2m-by-2m, that define the two dual operations. 

Proof. First note that U is unitary and U = U T 
Since 

UU T = - (7 <g) 7 -iZ® Z'f = iZ®Z' 
and, by conjugation, UW = —iZ Cg> Z', we find 

x i (g> Y i = (z ® z') (x T ® y T ) (z ® z') 
= uu^ (x T ®y T ) uuK 

Therefore 

u ] (x a ® y«) u = u T (x T ® y T ) 77 = {u ] (x ® y) u) t . 

□ 

Remark 1.4. We will see later that this lemma really has proven the isomorphism 

M JV (H) ® M m (H) ^ M 2N+2m (R) 

where H is the algebra of quaternions, with is a real C*-algebra. In contrast, the next lemma 
is a tautology, but a good place to discuss a critical conventions. 

Lemma 1.5. For all X G M 2N (C) and Y G M 2m (C), 

x T ®y* = (X(g)Yf. 
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Proof. We adopt the convention on identifying matrices in Mj(C) <g) M^(C) with matrices 
in Mj + fc that makes the association 



T <g> 



a b 




aT 


bT ' 


c d 


cT 


dT 



when k = 2. Using the other convention can change the argument of the Pfaffian, which is 
the only interesting thing about the Pfaffian. (The Pfaffian's magnitude is the the square 
root of the magnitude of the determinant.) With this convention, 

I®Z Z' 

where the Z are of appropriate size, and so 

X T <g> = - (/ ® Z) (X T <g> Y T ) (I®Z) = -(I®Z)(X® Y) T (/ ® Z) 

which, under the chosen identification, says 

x T <g> r s = (x ® . 

□ 

2. BOTT PERIODICITY AND SOFT REPRESENTATIONS OF THE ZERO SPHERE 

2.1. GUE in all dimensions. To construct the topological invariants, we define an operator 

d+l 

(2.1) B(H 1 ,...,H d+1 )=J2 H r®lr, 

r=l 

where the 7 r are a set of anti-commuting Hermitian matrices: {7 n 75} = 0. In the GUE case, 
then we simply choose the 7 r to be a set of matrices of the minimal dimension to provide 
d + l different 7-matrices. That is, for d = 2, we can choose the 7 matrices to be the three 
different 2-by-2 Pauli spin matrices, while for d — 3,4, the 7 matrices need to be at least 
4- dimensional. In the GOE,GSE cases, we will choose the 7 matrices as described later to 
make certain symmetries of B(H\, ...,H d+ x) more apparent. 

Assuming the matrices H r are a soft representation of the sphere, as in (jl.3|1.4p . then B 
is a soft representation of the sphere S°. This simply means that 

(2.2) B 2 w I and B ] = B. 

Thus, the eigenvalues of B are close to plus or minus one. 

The integer invariant we consider in the GUE case is one-half the difference between the 
number of positive and negative eigenvalues of this matrix, which for consistency with the 
GSE and GOE cases we give a second name, 

B(Hi, H d+ i) = B(Hi, H d+ i). 

When there are no zero eigenvalues we call this quantity bott(H 1 , H d+1 ). This index, 
called the Bott index, is a topological invariant, in that it does not change along any path of 
matrices H r which form a ^-representation of the sphere for sufficiently small 5 (the value of 
5 required depends upon dimension; for d = 2 we need 5 < 1/4). This is proven by showing 
that the only way for the Bott index to change is for an eigenvalue of B = B(Hi, H d +i) 
to become equal to zero. If this happens, then \\B 2 — I|| = 1. However, if a set of H r 
form a ^-representation of the sphere for sufficiently small 5, then \\B 2 — I|| < 1, giving a 
contradiction. One can similarly show that given two different tuples of matrices H r and K r 
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which both form ^-representations of the sphere, then if \\H r — K r \\ is sufficiently small then 
bott(i?i, .., H d+ i) = bott(i^i, .., K d+ i) (generalizing lemma 3.5 of [1]) by considering a linear 
path (1 -t)H r + tK r for t e [0,1]: 

Lemma 2.1. Suppose (Hi, H d+ i) and (Ki, K^i) are tuples of self- dual, Hermitian 
n-by-n matrices and suppose (Hi, ...,H d+ i) is a 5 -representation of the sphere with 5 < l/n, 
where n is defined to be d(d + l)/2 + 1. If 

d+l 



then 



\\H r - K r \\ < Vl-nS 

r=l 

bott{Ki, K d+ i) = bott^x, H d+ i) 



Further, if 



Lemma 2.2. If Hi, .H d +i are exactly commuting and bott(f/i, H d+ i) is defined, then 
bott(Hi,...,H d+ i) =0. 

These last two lemmas establish that if the Bott index is nontrivial and 5 is sufficiently 
small then the given tuple of matrices H r is not close to an exactly commuting tuple. 

Note that the Bott index is always equal to zero in the case that d is odd, since the matrix 
I eg) (7i72-.-7d+i) then anti-commutes with B. 

2.2. 2D GSE. We now consider the case of problems with time-reversal symmetry, where 
the Hamiltonian is in the GSE universality class. In this case, the presence of time-reversal 
implies that the Bott index above defined in the GUE class is trivial. However, there is 
another non-trivial index. This index represents an obstruction to finding localized Wannier 
functions which respect time-reversal symmetry. 

The procedure will be to show that in two dimensions, we can construct a unitary trans- 
formation that makes B(H\, H 2 , H 3 ) an anti-symmetric matrix. Then, the invariant that we 
consider is the Pfaffian of this matrix. The procedure in three dimensions is to construct a 
unitary transformation that makes B(Hi, H 2 , H 3 , H±) a real symmetric chiral matrix. Such 
a matrix is of the form 

' 

A T 



(2.3) 



If B 2 f« /, then A is approximately orthogonal. The invariant that we consider in this case 
is the determinant of this matrix. We notice a general pattern here: starting with a certain 
number of matrices (3 or 4) in a given symmetry class (GSE) we construct a single matrix 
in another symmetry class (anti-symmetric or chiral real symmetric, respectively). See table 
4 in [7], where the symmetry classes are arranged in a sequence given by Bott periodicity. 

We begin with the two dimensional case, reviewing the construction of [I], used in [9]. 
From a physical point of view, the existence of a unitary transformation that makes B 
anti-symmetric by is not surprising: the self-dual operation can be regarded as a time- 
reversal symmetry operation, and a similar time-reversal symmetry can be applied to the 7 
matrices used to construct B{H\, H 2 , H 3 ). We choose a time-reversal symmetry operation 
that makes those 7 matrices odd under time-reversal. Then, under these combined time 
reversal symmetries, B(H t , H 2 , H 3 ) changes sign; however, since there are two spin-l/2s, 
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the time reversal symmetry operator squares to unity and hence, up to a basis change, is 
equivalent to transposition. This is simply lemma [L3l above, as we will see. 

Definition 2.3. Let H r be self-dual and Hermitian. Define the matrix B(H±, H 2 , H 3 ) by 

(2.4) B(H h H 2 , H 3 ) = rtB(H x , H 2 , H 3 )U, 
where the unitary U is defined by 

(2.5) U=-={I + Z®a 2 ). 

v2 

This matrix B will be anti-symmetric, and so i times a real matrix, by lemma ll.3[ since 
Hi = H r while a\ = — oy. Still following PQ, we now define the index by taking the Pfaffian 
of this matrix B. Since B is Hermitian and pure imaginary so has real eigenvalues that occur 
in pairs, symmetric across zero. Thus its determinant is nonnegative and its Pfaffian is real. 

Definition 2.4. We define the index bott(ifi, H 2 , H 3 ) for self-dual matrices H r by 

(2.6) tott(#i, H 2 , H 3 ) = sgn(Pf (B(H U H 2 , H 3 ))), 

where Pf is the Pfaffian and sgn(x) = 1 for x > and sgn(x) = —1 for x < 0. If 
Pf(B(H u H 2 , H 3 )) = 0, the index bott(tfi, H 2 , H 3 ) is not defined. 

In pQ, it is also shown that this index is a topological invariant, in analogy to the GUE 
case: 

Lemma 2.5. Consider any continuous path of self- dual matrices, H r (s), where s is a real 
number, < s < 1. Suppose that for all s, the matrix B(Hi, H 2 , H 3 ) has non-vanishing 
determinant. Then, batt(i?i(0), H 2 (0), H 3 (Q)) = bott(#i(l), H 2 (l), H 3 (l)). 

Proof. The determinant of B(Hi, H 2 , H 3 ) is equal to Pf (B(Hi, H 2 , H 3 )) 2 . As long as the 
determinant does not vanish, the Pfaffian does not vanish and hence cannot change sign. □ 

and 

Lemma 2.6. If Hi,H 2 ,H 3 are self-dual and exactly commuting and bott(ifi, H 2 , H 3 ) is 
defined, then bott(f/i, H 2 , H 3 ) = 1. 

and 

Lemma 2.7. Suppose (Hi, H 2 , H 3 ) and (Ki,K 2 ,K 3 ) are triples of self- dual, Hermitian n- 
by-n matrices and suppose (Hi, H 2 , H 3 ) is a 5 -representation of the sphere with 5 < 1/4. 
If 

\\Hi - Ki\\ + \\H 2 - K 2 \\ + ||fl- 3 - K 3 \\ < VT^45 

then 

brtt(K u K 2 , K 3 ) = boTt(fT a , H 2 , H 3 ). 
This last lemma justifies calling this index an invariant. 

Given a soft-representation of the two-sphere with arbitrary complex matrices we can 
double these with their transposes to get self-dual soft-representation. The resulting Z 2 
index is determined by the parity of the original A index. 
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Theorem 2.8. For Hi, H 2 , H 3 a soft representation of the two-sphere, in Mjy(C), 



" Hi 







" H 2 







' H 3 








#i _ 


1 





H~ 2 _ 







H~ 3 



is a self-dual soft representation of the two-sphere and 

bott 
Proof. Recall 



" Hi 







' H 2 







' H 3 















H~ 2 _ 


1 





H~3 _ 



^_^Bott(H 1 ,H 2 ,H 3 ) 



/ o o -i/ 
/ il 
OHIO 
-il I 



and the sign of Pf (WBU) gives the Pfaffian-Bott index, where 



B = B 



Hi 
Hj 



H 2 
HJ 



H 3 
HJ 



Let 



Q 



l 



I 





-I 








i/ 





il 





7 





-I 


-il 





-il 






Notice 





'I -I' 




' il il ' 


(>/2 


il il 


I -I 



so Pf (WBU) = Pf (Q T WBUQ) . Since 



UQ 



and 



g T c/ f 









-I 








il 

















-I 


-il 











' I 

















I 














-il 





il 
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we find 



Q T U*BUQ = Q T U^ 



Hi 



H 3 
H x + iH 2 






HI 
H? - 1 HJ 



H 3 





Hi-i 





HI 





+ iH 2 





-H 3 





H? + iHl 







Hi 


-iH 2 







-H 3 





hi - mi 



-HI 



UQ 



-iHj + Hi 



iHl 






and so 



Pf (U j BU) = det(-f) det 



#3 

Hi + iH 2 






-H 3 
-Hi - %H 2 





Hi - %H 2 
-H 3 



UQ 






iHl 
iHl + Hi 

-Hi + iH 2 
H 3 





det 3 {Hi, H 2 ,H 3 



If bott(ifi, H 2 , H 3 ) = m then the spectrum of the 2N-by-2N matrix B (Hi, H 2 , H 3 ) will have 
2N — m negative eigenvalues, and so the sign of Pf (WBU) will be negative exactly when 
m is odd. □ 

2.3. 3D/4D GSE and 4D/6D/7D GOE. We now turn to the three-dimensional case. In 
this case, we choose a particular representation of the 7 matrices as 

(2.7) 71 = I® a* 



7i 

72 
73 



a, 



o,, 



a 



a 



a,,. 



74 = CT z vyvy, 

where a r are the Pauli spin matrices. This gives 4-by-4 7 matrices. Writing B = H r (g>j r , 
we have 4 different matrices textured together: the orbital and spin degrees of freedom of 
H r , and the two two-dimensional spaces used to define 7 r . 

The existence of a transformation making B chiral and real is not so surprising from the 
following physical point of view, again using lemma TOl and, similar to the two dimensional 
case, using the trick of picking a set of 7 matrices with appropriate behavior under an 
appropriate duality operation. The H r are self-dual, so that Hj = —ZH r Z. Now, consider 
the matrix Z' = ia y ® I. The j r are self-dual using Z' to define the duality: 7^ = —Z'^ r Z'. 
The matrix B is then a sum of tensor products of two self-dual matrices, H r and 7 r . So, 
B should be symmetric under a time reversal operation. However, given two spin-l/2s (one 
spin- 1/2 corresponding to the spin degree of freedom of the H r and the other being the first 
of the two sigma matrices used to define the 7 r ), we form a system with spin or spin 1, so 
up to a unitary the time-reversal operation is the same as transposition. This means that 
unitarily conjugating B gives us a symmetric matrix, and symmetric plus Hermitian means 
real symmetric. Next will then see why it is chiral. 

We see the change of basis to make this in in lemma 11.31 By that lemma, 

tfB(Hi,H 2 ,H 3 ,H 4 )U, 
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where the unitary U is defined by 

(2.8) U = ^={1 -iZ ® Z'). 

v2 

is real symmetric. 

We now show that this is chiral in an appropriate basis. First 

Lemma 2.9. The matrix WB(Hi, H 2 , H 3 , H 4 )U anti- commutes with the matrix I®I®I®a z , 
where the first two matrices refer to the orbital and spin degrees of freedom used to define 
H T , and the last two refer to the two two-dimensional space used to define j r . 

Proof. The matrix / <8> I (g) / <8> o~ z anti-commutes with U, so we must show that B anti- 
commutes with I ®I <S)I ®a z . However, this follows since I ®I ®I ®o~ z anti-commutes with 
all the matrices j r . □ 

Thus, if we write W B{H\, H 2 , H 3 , H 4 )U as a block matrix, with the two blocks corre- 
sponding to the positive and negative eigenvalues oiI®I®I®o~ z , we have a symmetric 
chiral real matrix. 

For numerical purposes, of course, it is convenient just to compute the upper-right-hand 
block of Eq. (12. 3p . since the lower-left hand block is related by transposition. This block can 

be expressed as 

-(/ - i<Ty)H X (I + io v ) + %B{H 2 , F 3 , H 4 ), 

and it is essentially this block that we take as B in this case. To be consistent with what we 
coded, we define this as follows. 



Definition 2.10. Let H r be four self-dual and Hermitian matrices. Define the matrix 

(2.9) B(H U H 2 , H 3 , H 3 ) = ^ H r <g> v)j U, 

where the unitary U is defined by 

(2.10) U= - 1 ={I + Z®a 2 ). 

v2 

and 

(2.11) v x = I, v 2 = icr x , u 3 = ia y , u 4 = ia x . 
We now define the index for this problem as: 

Definition 2.11. We define the index bott(ifi, H 2 , H 3 , H 4 ) for self-dual matrices H r by 

(2.12) boTt(i?x, H 2 , H 3 , H 4 ) = sgn(det( J B(F 1 , H 2 , H 3 , H 4 ))). 

We summarize in Table [2] our method to deal with approximate representations of S d by 
self-dual matrices, meaning Hj. — H r , HI — H r and 

d+l 
r=l 

in dimensions 2, 3 and 4. In all cases, 

u = -L(i + z®z') 
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Relations on 

v \i ■ ■ ■ i v d+l 


Definition of B = B(H U . . . , H d+l ) 
and choices for v\, . . . , Vj+i 


Relations 
on B 


Scalar 
valued 
function 


all d 


4 = k 1 

v\v s = —v\v r 
for r s. 


/d+i \ 

B = H r ® 


B^B « / 




d = 2 


t 

K = V r 




7~~) f 7~> 

B ] = B 


2 signature 


d = 4 


I/T = V r 


£7 X (J y 

0^ j ' [ -Gy _ 
— IGy 
IGy ' C 


r o i 

' [ G Z \ ' 

3 (Xj, 1 


B ] = B 


^Signature 


d = 6 


V\ = Vr 




B^ = B 


| Signature 


d = 8 


v\ = V r 




B^ = B 


^Signature 



Table 1. GUE approximate representations of S . The signature of B leads 



to an integer invariant. 



where 



Z' 



I 
-I 



in the appropriate size. Table |3] deals with GOE case. Notice that in some dimensions the 
Bott matrix B has different symmetries than the original matrix, so real H r can lead to self- 
dual B. Tables [5] and [6] show how these relations lead to i^-theory elements for an algebra 
of matrices over the reals of the quaternions. 

Table [1] shows the way to deal with approximate representations of S d in the GUE case, in 
dimension 2 and higher even dimensions. Table H] shows how these relations lead to i^-theory 
elements. 

The following theorem tells us that we a dealing with a real, or purely imaginary, invertible 
matrix with the correct symmetry so that we can compute a real Pfaffian, a real determinant, 
or perform a positive/negative eigenvalue count of the real eigenvalues. 

Theorem 2.12. Suppose Hi, ... , Hd+\ are 2N-by-2N matrices with Hj. = H\ = H r , 

d+l 

^H r -I <5 

r=l 

and \\[H r , H s ]\\ < 5 for r ^ s and some 6 > 0. Suppose v\, . . . , i>d+i are matrices in M 2m (C) 
so that 



v. 



v\v s = -v\v r (r ± s) 
for r s. Let Z^ be the matrices defining the dual operations on the M 2 fc(C) and define 



U 



V2 



(I + Z 



N 
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Relations on 

Z/i,...,Z/ r 


Definition of B(H\, . . . , H d+ i) and 
choices for ui, . . . , v d +\ 


Relations 
on £ 


Scalar 
valued 
function 


all a 


4 = ^ 
for r ^ s. 


/d+i \ 

B r ® v r 1 (7 






d = 2 


v\ = v r 
v\ = -v r 




B^ = B 
B T = -B 


Pfafnan 


a = 3 


& t 
z/ H = — v ] 

r r 


I,io x ,io y ,io z 


Z5 1 = ±> 


Determinant 


d = 4 


V), = Vr 
v\ = V r 


[ cr x I Oy 

<7 X j ' [ -Oy _ 

— io y 

%Oy ' c 


' 

' [ a z y 

3 (T^ 

-y 


B ] = B 
B T = B 


\ Signature 


d = 8 


v\ = V r 
T 

V r = V r 




B* = B 
B^ = B 


^Signature 



Table 2. GSE approximate representations of S d . The signature of B leads 
directly to an integer invariant, while we must take the sign of the determinant 
and Pfafnan to get ±1 and so a Z 2 invariant. 





Relations on 
v u ...,v r 


Definition of B(H\ ) . . . , H d+ i) and 
choices for ui, . . . , v d+1 


Relations 
on B 


Scalar 
valued 
function 


all d 


4 = 

v\v s = -v\v T 
for r^s. 


d+l 

B = H r ® v r 

r=l 


B^B « 7 




d = 4 


v\ = v r 
v\ = v T 


l <J X 1" Oy 

a x \ ' [ -Oy 
— io y 
ioy ' c 


r <?z ■ 

' [ z 

3 Oy 
-y 


B^ = B 
B* = B 


^Signature 


d = 6 


v\ = v r 

T 

V,. = -V r 




= B 
B T = -B 


Pfaflian 


d=7 


u T = -i/t 

r r 




B^ = B T 


Determinant 


d = 8 


v\ = v r 

T 

V^ = V r 




B^ = B 
B T = B 


\ Signature 



Table 3. GOE approximate representations of S d . The signature of B leads 
directly to an integer invariant, while we must take the sign of the determinant 
and Pfaflian to get ±1 and so a Z 2 invariant. 



'd+l 



B = B(H l7 . . .,H d+1 ) = \J2 H r ® "v ]U. 



K r=l 
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(1) We have 

\\B^B-l\\ < (d + 2)5. 

(2) // i/t = v r for all r then B ] = B. 

(3) // v\ = u r for all r then B T = B. 

(4) If v \ = -u r for all r then B T = —B. 

Proof. (1) We use the fact that U is a unitary and we find 

UB*Brt = J2 H r®I + J2 H »\ ® ^ z 



r<s 



SO 



and 



UB ] BU ] -1= f H r - 1 ) ® / + ^ [#r, #.] ® ^ 



Ists-ill < 



+ J2\UHr,H,]\\\\4 



r<s 



The v r are unitary, so = 1 and we have the desired upper bound. 

(2) From the definition of B we quickly obtain 

'd+l 



so i>l = v implies B^ = B. 



(3) and (4) Lemma [1.31 implies 



,r=l 



'd+l 



B T = £/t [J2 H r <g> z/M U 



,r=l 



□ 



and so v\ = ±v r implies B T = ±B. 
Theorem 2.13. Suppose Hi, ... , H d+ i are n-by-n real symmetric matrices with 

d+l 
r=l 

and || [H r , H s ]\\ < 5 for r ^ s and some 5 > 0. Suppose u±, . . . , Vd+\ are matrices in M 2m (C) 
so that 



v\v s = -v\v r (r ^ s) 



for r s. Let 



d+l 



B = B(H 1 , . . . , H d+ i) = 22 Hr 



r=l 



(1) We have \\B*B - l\\ < (d + 2)5. 

(2) // v\ = v T for all r then £t = B. 

(3) // U J = v r for all r then B* = B. 
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TV" / A \ 


For K n (C). 






Restrictions 
as given 


Equivalent 
conditions 


Numerical invariant 


n = 
0,2,... 


xt = X 


Hermitian 


|Sig(X) in Z 



Table 4. K n (A) built using invertible elements in M n (A) for a C*-algebra 
A. Here Sig(X) denotes the inverible matrix X, the difference between the 
number of its positive eigenvalues and the number of its negative eigenvalues. 
The rows repeat in the pattern as Bott Peridicity has order two in the complex 
case. 



K n (ft(A)) 


For A n (R) = 


= A n (K(C,id)). 






Restrictions 
as given 


Equivalent 
conditions 


Numerical invariant 


n = 


X T = X 

xt = X 


real and 
symmetric 


|Sig(X) in Z 


n = 1 


x T = xt 


real 
orthogonal 


sgn (det(X)) in Z 2 


n = 2 


X T = -X 

xt = X 


pure- 
imaginary 
and hermitian 


sgn(Pf(X)) in Z 2 



Table 5. The first three A-theory groups Mfc(R) of the are summarized. 
Here A refers to a generic real C*-algebra. 



(4) If vj = -v r for all r then = -B. 

Proof. The proof is nearly the same as before, but we now use lemma 11.51 □ 

In addition to the Pfaffian, we use the determinant and signature. 

Definition 2.14. If X is an invertible, Hermitian matrix, its signature Sig(X) is the number 
of positive eigenvalues of X minus the number of negative eigenvalues of X. 

We will see below that the sign of the Pfaffian is identifying a class in K2(M), the sign of 
the determinant is identifying a class in ^(R), while the signature is identifying a class in 
A"o(M). The reason this A-group appear is that the interesting i^-theory of the <i-sphere ends 
up in fQ-d which results in an invariant in i^ 8 _^(H). Rather than attempt to understand 
directly K 6 , K 5 and K± of the quaternions H, we use isomorphisms 

K n (U) S A n _ 4 (H <8) H) = K n ^(R) 

which means that in the GSE case, a system in dimension d has an invariant in Kt-dQH). As 
the A-theory of R is 

Z,Z 2 ,Z 2 ,0,Z, 0,0,0 

in degrees 0, 1, . . . we are getting, as expected, invariants in Z 2 in dimensions 2 and 3, and 
an invariant in Z in dimension 4. 
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K n {U{A)) 



n 



n 



n 



K n {W) = K n (5ft(M 2 (C), 



Restrictions 
as given 



X« = -X 

xt = X 



x* = x 



Equivalent 
conditions 



pure- 
imaginary 
and hermit ian 



real 
orthogonal 



real and 
symmetric 



Numerical invariant 



iSig(X) in Z 







Table 6. The first three if-theory groups of the the quaternions Mjy(H) 
(realized as matrices respecting the dual operations) are summarized. The \ 
is needed instead of | to account for the doubling of eigenvalues in self-dual, 
Hermition matrices. 



3. Mathematical Problems 

There are several interesting mathematical problems which arise from this construction. 
This section will be useful even if one is only interested in numerical applications of this 
techniques, since an understanding of what has been proved and which directions of the 
implications are known will be very useful. The first five problems are purely problems in C*- 
algebra, having to do with matrices, while the last problem involves fermionic Hamiltonians. 

In some of these problems, we consider the question of stable limits for matrices. This 
is very similar to the idea of the stable limit for free fermion systems. The approximation 
problem in the stable limit, given matrices H r which are a soft representation of S d , is 
whether one can find exactly commuting matrices J r , with J2 r J? = suc ^ that the set of 
matrices {H r © J r } can be approximated by a set of exactly commuting matrices H' r . 

We begin with 

Problem 3.1. Consider d + 1 matrices H r , belonging to one of the 10 different random 
matrix theory universality classesj38], giving a soft sphere S d . Determine the topological 
invariants in the stable limit. 

We conjecture that 

Conjecture 3.2. If the H r belong to classes GUE,GOE, or GSE, then the solution to the 
above problem is given by the same results as in [HI [7] in the free fermion case. See, for 
example, table 4 in [7J, where the dimension along the top is the dimension d of the sphere, 
and the universality class along the left is the universality class of the matrices H r . The 
precise conjecture is that for soft representations of the sphere in dimensions 2 through 9, 
the only topological invariants in the stable limit are invariants described in Tables (H |2] and 

El 

In fact, as we stress, only in some cases is this conjecture what we need to study matrix 
invariants of free fermion systems. Given free fermion systems in the GOE,GUE,GSE classes, 
by forming projected position matrices H r we obtain matrices H r which are in the same 
universality class as the original Hamiltonian. In the previous section, we explained how to 
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construct invariants in two particular cases, the GSE case in d = 2 and d = 3. Later, we 
consider this from a more general K-theoretic point of view, giving a calculation of invariants 
of these three classes in all dimensions. The conjecture above is essentially equivalent to the 
conjecture that the invariants we have computed comprise all of the invariants in these classes 
in the stable limit. 

Problem 3.3. Show that a nontrivial value of the invariant computed from the B(Hi, i?<z+i) 
(either a nontrivial integer or Z 2 invariant) implies that the matrices Hi,..., Hj+i cannot be 
approximated by exactly commuting Hermitian matrices of the given symmetry class, GUE, 
GOE, or GSE. 

This result has long been known in the case of d = 2 without symmetries, and we have 
given a proofpQ in the time reversal invariant case in d — 2. The proof in other cases 
(other symmetry classes or dimensions) is very simple and completely analogous to the proof 
in the other known cases: any other set of matrices H' r with \\H' r — H r \\ small gives rise 
to an B(H[, H' d+1 ) with B(H[, H' d+1 ) — B(Hi, Hd+i) small. However, the matrix 
B(Hi, ...,Hd+i) is close to a projector, so its eigenvalues are far from zero, so invariants of 
the matrix such as number of positive eigenvalues, determinant, or Pfaffian, do not change 
under small changes in B. So, this problem is solved. 

Problem 3.4. Show that a trivial value of the invariant computed from B[Hi, H^+i) 
implies that the matrices can be approximated by exactly commuting Hermitian matrices of 
the given symmetry class, either GUE, GOE, or GSE, possibly in the stable limit. 

This problem has been solved in the case of Hermitian matrices with no other symmetries. 
In the case of d = 2, this is solved in [15] and the solution is reviewed in [TJ. In the case 
of d = 2, this result holds without going to the stable limit; that is, it is not necessary to 
add on trivial degrees of freedom. This relies heavily on the fact that two almost commuting 
Hermitian matrices can be approximated by two exactly commuting Hermitian matrices, 
known as Lin's theorem [TBI E]- m more than two dimensions, this has been solved in the 
stable limit [12] . 

Problem 3.5. Show that for each dimension and each symmetry class and for any given 
choice of the invariant in that symmetry class, we can construct, for any sufficiently large 
integer N, d + 1 different N-by-N matrices H r in the appropriate symmetry class, with 
|| [H r , H s ] || converging to zero as N tends to infinity and || 'Yl ia H^ — 7|| also converging to 
zero as iV tends to infinity and such that, for any N, the matrices have the given value of 
the invariant. 

In fact, in the case that the invariant is an integer, it suffices to find matrices for which 
the invariant is equal to +1 and to —1, as any larger integer invariant can be obtained 
by the direct sum of matrices with invariant +1 or —1. This problem has been solved in 
the case of d — 2 with time-reversal symmetry or with arbitrary Hermitian matricespQ. A 
general solution would be obtained if one could solve this problem for every dimension. One 
route to obtaining this mathematical result is, of course, to construct the desired matrices 
by constructing an appropriate free fermion system on a sphere and then constructing the 
band projected position matrices following our procedure. This would require solving the 
next problem, that the matrix invariants agree with the fermionic invariants: 

Problem 3.6. Show that the matrix invariants we compute agree with the classification of 
free fermions of |8j, so that for each fermi Hamiltonian we obtain the same invariant. 
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We now sketch a possible approach to showing this in the Z case (the Z 2 case would 
be similar, so we just sketch one case). We would begin by showing that the invariants 
agree for any single system H + \ with the invariant equal to +1 and any other system if_i 
with the invariant equal to —1. We use the term "matrix invariant" to denote the invariant 
we compute and we use the term "free fermion invariant" to denote the invariants from the 
classification of |8J. Once this is done, we can show that the invariants agree for all systems as 
follows: let H be a given free fermi Hamiltonian. Suppose H has a free fermion invariant 
equal to m. Suppose m > 0, without loss of generality. Then, consider the Hamiltonian 
H = Hq®H-i®...®H^i, where we add m copies of H-\. Then, by definition of the invariants 
of [8], this Hamiltonian can be continued to a trivial Hamiltonian, so this Hamiltonian H 
has localized Wannier function. Thus, H has a matrix invariant equal to 0. However, the 
matrix invariants we study are additive under direct sum. That is, let B({H r }) denote the 
index computed from a given set of matrices H r . Given any matrices H r and H' r , we have 

(3.1) B({H r }) + B({H' r }) = B({H r © H' r }). 

Choose H' r to be the band projected matrices obtained from the Hamiltonian if_i © .. © if_i 
and let H r be the band projected matrices obtained from the Hamiltonian H . Then, the 
right-hand side of this equation is equal to zero since H has matrix invariant equal to zero, so 
Eq. (13.11) implies that B({H r }) = m. So, it suffices to show equality in only one non-trivial 
case. 

4. ID Systems and the Polar of a matrix 

4.1. Soft one-torus and one-sphere representations. In one- dimensional GOE, GUE 
or GSE systems (non-chiral) systems, there are no obstructions to gapped and local Hamil- 
tonians being deformed to trivial Hamiltonians. We can use a basic matrix function to fix 
the one band compressed periodic observable to make it an actual unitary, diagonalize that 
unitary and so find local Wannier functions. We go through this simple case in detail as we 
develop machinery used in the higher dimensional cases. 

The simpler nature of ID systems is related to the fact that in C*-algebras and matrix 
theory, the relations associated to a one-dimensional space are "stable" in a way that fails 
for the relations associated to higher-dimensional spaces. For example, consider a matrix X 
such that 

(4.1) ||X|| < 1, X ] 

It is trival to see that with Y — A [X^ + Xj we have Y ~ X and 

(4.2) ||y|| < i, yt = y. 

This same trick works in any C*-algebra. As the eigenvalues of Y are real and of magnitude 
at more 1, we think of equation ( 14.4f) as being associated to a line segment. It is so easy 
to fix "approximate Hermitian" operators to be Hermitian that generally one only studies 
inexact relations between Hermitian matrices, or relations between general matrices. 

This approximation of almost Hermitian X by exactly Hermitian Y is an example of a 
class of problems called "stable relations". We have a set of matrices that almost obey certain 
constraints, such as a matrix being almost Hermitian, and we ask whether the matrices can be 
approximated by matrices that exactly obey the given constraints. Our discussion in most 
of this paper is centered on the question of stable relations involving almost commuting 
matrices: giving a set of almost commuting matrices, can they be approximated by exactly 
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commuting matrices? However, in parallel with this question of almost commuting matrices, 
we often have to deal with problems such as replacing almost unitary matrices by exactly 
unitary matrices, so it is worth mentioning the more general problem of stable relations. 

Moving to the square, it is a difficult theorem [13] that given a matrices X 1 and X 2 such 
that 

(4.3) ||X r || < 1, X\ = X r , X X X 2 « X 2 X X 
there will exist Y r matrices with Y r « X r and 

(4.4) ||Y r || < 1, Y} = Y r , Y X Y 2 = Y 2 Y X . 

Properly stated, with uniform norm conditions, this is false in general C*-algebras. 

Some relations are not obviously associated with any topological space, but give interesting 
stability results, even for matrices. For example, given a matrix X with 

||X|| < 1, X n « 

there is a matrix Y with Y ~ X and 

||Y|| < 1, Y n = 0. 

This is true, but not for easy reasons, in general C*-algebras [16] . 

What we need for the study of ID systems is a way to find for matrices Xi and X 2 with 

(4.5) Xl = X r , X X X 2 « X 2 X U X\ + X 2 2 ^I 
matrices Y± and Y 2 with Y r Y r and 

(4.6) Y} = Y r , Y X Y 2 = Y 2 Y h Y? + Y 2 2 = I. 

As the scalars that satisfy equation (14.61) are just sin(0) and cos(6>) we consider equation 
(14.61) as being associated to the circle. Such a ID "stable relations" problem has an easy 
solution. We will verify that this solution respects the needed symmetries so that when we 
fix the one band compressed periodic observable and so find local Wannier functions, those 
will have the correct symmetry as well. 

The only geometry to consider in ID is the one-circle/one-torus. There is a distinction in 
the equations used to describe the same space that becomes important when we compress 
the position operators. We consider a lattice on a torus of radius L and we use 6 to denote 
angle on the torus. We define 6{i) to be the angle 6 of site i, and let n be the number of sites. 
In analogy to the two-sphere case, we define to be a diagonal matrix, with On = 6(i), in 
M n (C ?l ). We define band projected position matrices with 

Pexp(z&)P = Q^ U Q °^JQ\ 

where Q is a unitary (suppressed above) for the change of basis that puts U in a canonical 
block where it can be extracted numerically. As before, P almost commutes with exp(i@) 
so that U is almost unitary. We call this U a soft-representation of the one-torus, meaning 
only that 

H^tz-Jll «0. 

(In the infinite dimensional case we would need to add the relation \\UW — J < 1.) 
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If we were in the GSE case, we would be working with 2N sites, where ej and e N+ j occupy 
the same point but typically correspond to spin up and spin down. We would insist on H 
and being self-dual, so 6 would be diagonal with diagonal 

9(1),...,0(N),9(N + 1),...,9(2N). 

Observables that are Hermitian are more standard, so returning to the GUE case we 
introduce 

X 1 = Lcos(0), X 2 = Lsin(0). 
Then we can band-compress these and define H x and Hi by 

px 1 p = lq(1 °W 



H x 

PX 2 P = LQ [° n ° ) Qt 



H 2y 

These are what we call a soft-representation of the one-sphere which means 

H\ = H u H\ = H 2 , Hi + Hi « /, [H u H 2 \ « 0. 

Since 

G(S 2°H 1 ) Qj = 2PCOS{@)P = Q ^ U + Ui)&> 
and using a similar equations regarding H 2 , we have the expected equations 

ITx = ± (tf + E/+) , i/ 2 = £ (C/ - £/+) 

that allow for an easy translation between the two situations. 

We will obtain Wannier functions by a three step process: adjust Hi and H 2 so they exactly 
commute and exactly square-sum to the identity; jointly diagonalize these new matrices to 
produce a basis of common eigenvectors; use Q to move these common eigenvectors to the 
larger space where they become a basis of the low-energy subspace PC™. As they will be 
approximate eigenvectors of the original observables, they will be somewhat localized. 

If the Wannier functions are to have the correct symmetries in the GSE and GOE cases, 
then we must choose Q carefully and do the simultaneously diagonalizing in a way that 
respects the needed symmetries. We address this after explaining how soft representations 
are adjusted to be exact representations. 

Theorem 4.1. Suppose H 1 and H 2 are matrices such that HI — Hi, H\ — H 2 , 

11^ + ^1-/11 < 5 

and 

\\[H U H 2 ]\\<S. 

If 8 < 0.6 then there are matrices K 1 and K 2 so that K* = K r and 

\\K r - H r \\ < 25 

for r = 1,2, 

K\ + K\ - I < 5 

and 

[K U K 2 ]<5. 

If the matrices H r are self-dual, then the matrices K r may be chosen to be self-dual. If the 
matrices H r are real, then the matrices K r may be chosen to be real. 



TOPOLOGICAL INSULATORS AND C*-ALGEBRAS: THEORY AND NUMERICAL PRACTICE 25 



polar (diag (Ai, . . . , A n )) = diag 



The proof, given in full below, works with the approximate unitary U = Hi + iH 2 . This 
seems an odd choice in the case of real matrices, but there are enough symmetries at hand 
to be sure we end with real matrices in that case. An approximate unitary can be dealt with 
using what is perhaps the most important of all matrix functions (the idea of a function 
of a matrix is sometimes referred to as "functional calculus"), the mapping of an invertible 
matrix to its polar part. 

4.2. The polar part of an invertible matrix. It is well known that for X an invertible 
matrix, there is a unique way to express it as a product X = UP with U unitary and P 
positive (meaning what applied mathematicians call positive semidefmite), and if X is real 
we find P is real and U is real orthogonal. Lesser known is the fact that if X is self-dual, 
then we find U is self-dual. A formula for U is U = X(X'X)~3. 

Definition 4.2. Given an invertible matrix X, the polar part of X is 

polar(X) = X(X*X)-*. 

A quick explanation for the name is 

A. A 

|Ai|"'"'|Ai| 

We can calculate polar (X) several ways. In the implementation section we discuss how 
Newton's method can be used to quickly compute the polar part of a matrix. Another 
method is to diagonalize X, so X = WDW\ and use the formula 

po\ar(WXW^) = Wpo\ar(X)W^ 

which holds for all unitaries, and the apply the scalar function A h- > A/|A| to the diagonal 
elements of D. 

Lemma 4.3. Suppose X is an invertible matrix in M n (C) with 

\\X ] X-l\\ < 5. 

I/O < 5 < 0.6 then 

||polar(X) -X\\ < 5. 
Proof. If X is invertible and has polar decomposition X = UP then 

||polar(X) -X\\= I - (X ] X)~ 

= max j|A~5 -l| | A e a (X+X) } 

< max ((1 - 5)^ - 1, 1 - (1 + S)~ 

For the range of 5 under consideration, this quantity is less that 5. □ 

Lemma 4.4. If X is an invertible, self-dual matrix in M 2 jv(C) then polar (X) is a self-dual 
unitary. If X is an invertible matrix in M 2 7v(IR) then polar(X) is real orthogonal. 

Proof. All parts of this are standard except the claims about the symmetries. We prove 
these, and more, in Theorem 19.51 □ 



Now we prove Theorem 14.11 
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Proof. Suppose H x and H 2 are self-adjoint, 

\\Hl + H 2 2 -I\\ < 5 

and || [H h H 2 ]\\ < 5. Let X = H l + iH 2 so that H x = \ (X* + X) and H 2 = -f (X* - X) . 
Let 

[/ = polar (X) 

and define ifi and K 2 by 
and 

K 2 = -| (17+ - IT) . 

These are evidently self-adjoint. Since U is unitary it commutes with W = U^ 1 and so K\ 
commutes with K 2 . Also 

K\ + Kl = \ (U ] U ] + 2U ] U + UU)-\ {U ] U ] - 2U ] U + UU) = I 

so we have an exact representation of the one-circle. As to the amount we have moved the 
H r to get the K r we estimate 

\\K r - H r \\ < - - X f || + - \\U - X\\ = \\U - X\\ . 

Lemma [4.31 gives us the estimate 

\\U-X\\ < \\X f X-l\\ 

= \\H* + HI-I + iH x H 2 - iH 2 H x \\ 
< 25 

and so \\K r — H r \\ < 25. 

If H\ = Hi and H\ = H x then A« = X. By lemma S3 = U and so K\ = K r . If H x 
and H 2 are real, then as they are self-adjoint they are symmetric. By lemma I4.4[ U T = U 
and so Kj — H r . As the K r are Hermitian, they are real. □ 

4.3. A structured spectral theorem. There are various versions of the spectral theorem, 
involving commuting matrices, that or Hermitian or at least normal, real or complex. We 
need a version that is not well known, involving self-dual matrices that are self-adjoint. 

Lemma 4.5. If X in M 2 jy(C) is normal and Xv = Av then 

X» (TV) = A (TV) . 
If X is any matrix in M 2 at(C) and Xv = Av then 

(Tv) t X» = A(rv) t . 

Proof. The ordinary spectral theorem tells us = Av. Conjugating this we discover 

A T v = Av. Since — ZX^Z = X T this means — ZX^Zv = Av which is equivalent to X^Zv = 
\Z\. Recalling from ( 11 . 14}) the definition of T, we finish by multiplying by —1. 

For the second claim, we start with the transpose of the eigenequation, v T A T = Av T . 
This implies -v T ZX«2 = Av T which solves to (v T Z) X i = A (v T Z) . □ 

The following is the self-dual finite-dimensional version of the spectral theorem. We assume 
familiarity with the complex and real versions. 
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Theorem 4.6. If Hi, . . . , H^ are commuting self-dual, Hermitian matrices in M 2 jv(C) there 
is a symplectic matrix U so that U^H T U is diagonal for all r, where the diagonal matrices 
are of the form 

' A r 
A r 

Proof. A finite set of commuting matrices will has a common eigenvector, so let v be a unit 
vector so that 

H r \ = A r v 

for all r. By lemma 14.51 we conclude 

H r (TV) = A r (TV) 

as well. Choose any symplectic unitary U\ so that Uiei = v. By lemma 11.11 we also have 
Uie N+ i = TV. Let K r = U\H r U x . Then 

K r ei = \ r e± 

and 

K r e N+ i = Aj.ejv+1- 
As the K r are self-adjoint, we conclude that in terms of N-by-N blocks, 



Unless N = 1, and we are done, we form 

B r D r 

These form a self-dual commuting family of Hermitian matrices. As simple induction now 
finishes the proof. □ 

4.4. Structured band-compressed position operators. Recall we had commuting Her- 
mitian matrices X\ and X 2 representing the position observables, and we adjusted their 
band-compressed versions PX r P in two ways to define Hi and Hi- We rescaled them to 
account the physical size of the lattice, and we changed basis so as to remove blocks of zeros: 





A r 



























A j. 








B r 





D r 



PX r P = LQ 



H r 



Q 



Abstractly P is /("H) for H the Hamiltonian and / the indicator function for the set 
(— 00, Ep\. However, we describe this more concretely to facilitate the later discussion of 
the numerical method, and to clarify a subtle point that arises when selecting Q in the GSE 
case. We use the eigensolver described in Section [941 so that the self-duality of the Hamilton- 
ian is reflected in the diagonalization and the computed matrix P will also be self-dual. The 
numerical error that accumulates in other eigensolvers could destroy the expected self-duality 
of P if the Hamiltonians has multiple eigenvalues very close to each other. 

In the GUE case we can work with any spectral decomposition of the Hamiltonian H with 
eigenvalues 

Ai > A 2 > • • • > A„ 
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and associated eigenvectors 

qi, ,Qn 

We set Q = [qi, . . . , q n ] to obtain the unitary the diagonalizes W to diag(Ai, . . . , A n ). If I 
represents the identity of size no- by-no, where n — n is the largest index with A„_„ > Ep, 
then 

" 
/ 



P = Q 



Q ■ 



Therefore 



(which is a partial isometry) and 



PQ = Q 

















-Q^PX r PQ 



1 

L 












QX r Q 








That is, H r is the bottom-right rio-by-rio block in j^QX r Q. 

We know there are K\ and K% close to H\ and if 2 that are exactly commuting and exactly 
square-sum to one. We discussed in pQ how this produces Wannier functions, localized in a 
somewhat weak sense. Here are illuminate the translation from the K r in M„ (C) to a basis 
for the low-energy subspace of C n . 

Let 6 be the larger of \\Hf + Pi\ — 7|| and ||[-ffi, -H2HI • That 5 will be small is discussed 
in the first section. By Theorem 14.11 we may assume K r is within 25 of H r . We can jointly 
diagonalize K\ and K2, with basis vi 

(4.7) K r v 

The corresponding basis of PC n is 

Q 



v„ and 



(4- 



Q 



Vl 







'no 



These are the desired Wannier functions. From ( 14. 7ft we derive 



and 



and finally 



PX r PQ 



Vl 





KQ 



Vl 





X r Q 



Vl 





Vl 





Precise estimates are possible, see [T], but the point is that an approximate eigenvector for 
the diagonal matrices X\ and X 2 , with approximate eigenvalues a and (3, must have small 
coefficients at sites far from (a, (3). 

In the GOE case we again can use a generic eigensolver to find a real orthogonal matrix 
Q that otherwise works as above. We obtain real Wannier functions, as desired. 

In the GSE we need a symplectic diagonalization of "H, meaning eigenvalues 

Ai ^ A2 " " ' A jy 



and their doubles 



A 



JV+l 



Ai,A 



N+2 



A 



2- 



A 



2N 



A 



N 
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and associated eigenvectors 

qi,...,qAr,Tqi,...,Tqjv 
Let Nq be the largest index with Aat_at > Ep. We assemble these eigenvectors to form a 
unitary Q as follows (in essence, we write first all the eigenvectors ^ with Aj > E F , then 
write the corresponding Tqi, then do the same for the eigenvectors qi with A« < Ep), 

Q = [qi> • • • , Qat-atq, Tqi, • • • , Tqjv_iv , qjv-iv +i, • • • , qv, Tqjv-jv +i, • • • , Tqv] . 
This matrix Q satisfies a symmetry similar to being symplectic, 

Z 

Z 



Q 



= ZQ 



where the bottom Z is of size 2N -by-2N . For matrices A and B, of appropriate sizes, 



Q 



A 
B 



Q 



Q 



Ai 
B» 







t 



The projection P will be self dual. We can even see this using Q to compute it: 

~ Ai 



H = Q 



A 



N 



Ai 



A 



N 



so 



P = Q 



and 



Thus also X 1 and X 2 are self-dual, and 




Hi 



Q 




H r 



= -(QQipX r PQQy 

Lj 



Q 



o o 

Hi 



which at last gives us ifj — H r . 

We can therefore find K r that are self-dual, and so obtain a common basis of the form 

Vi,...,Vjv ,Tvi,...,7vjv . 
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The Wannier functions are then 





Q 



Vl 



...,Q 







,Q 





TV! 



■ ,Q 







These have the desired structure, as 



T Q 





v, 



-ZQ 



Q 







Q 



Tv, 



5. Torus and Other Geometries 



We can consider lattice systems on other topologies such as the torus. In general, our 
procedure, discussed below is to map the matrices on the torus to matrices on a sphere and 
then compute the sphere invariants using the techniques described before. However, we begin 
with a particularly nice technique that is available in the GUE case on the two-torus. This 
technique will be useful later in relating the index to the Hall conductance. Unfortunately, 
we do not have access to such a simple formula for the index on the torus in the GSE case 
or in the GUE case on higher dimensional torus and so we have to resort to the mapping. 

We begin by describing this index on the torus. We then explain in the next subsection 
how to map the matrices from the torus to the sphere to define an index that way, and then 
theorem 15.31 shows that the two indices are the same in the GUE case. 

We consider a lattice on the torus, and we use 61,62 to denote angles on the torus. We 
define 6\{i) to be the angle 61 of site % and 6%{%) to be the angle 6*2(2). In analogy to the 
sphere case, we define G a for a = 1, 2 to be diagonal matrices with (O a )u = 6 a (i). We define 
band projected position matrices with 



(5.1) 



Pexp(i6i)P 
Pexp(iQ 2 )P 





Z7i 



u 2 



As before, P almost commutes with exp(zG a ), so that U\ and U 2 almost commute with each 
other and are both almost unitary. We call this a soft-representation of the torus. 
We now define the torus index. Let 

(5.2) tr(log(t/iZ7 2 t/Jt4) = r + 2mm, 

for some real numbers r, m. However, in fact m is an integer, because 

trflogf^i^C/M) = 

log(det(£A[/2t/M)) = log(det(|f/i| 2 )det(|Z7" 2 | 2 )), and det(|£/i| 2 )det(|£/ 2 | 2 ) is real. In taking 
the log in Eq. (15. 2p . we place the branch cut along the negative real axis. We define the 
index to be the integer m. 

We can join U r to polar(£7 r ) by a continuous path, and if 

||Etfl7 r -l|| <5, \\[U U U 2 ]\\ <5 

for sufficiently small S, then this torus index cannot jump to the next integer, and the 
commutator of the polar parts will stay small. This is important because the prior work 
on this index [17] was on almost commuting unitary matrices, while here we consider the 
problem of almost commuting matrices which are almost unitary. The point is that the index 
previously considered for almost commuting exactly unitary matrices also works for almost 
commuting almost unitary matrices. 
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5.1. Prom the soft torus to the soft sphere. We now describe various possible functions 
which can be used to map from the torus to the sphere, and how to use those to give an 
alternate way of computing the index for torus systems. Any continuous function 7 : T 2 — > S 2 
can, in theory, be used to determine a conversion 

(Ui,U 2 )^(Hi,H 2 ,H z ) 

that takes an approximate representation of the torus to an approximate representation of 
the sphere, so that we can then compute the index of the resulting representation of the 
sphere. 

A minor problem is that such a conversion is not uniquely determined by the function 7. 
More substantial issues have to do with i^-theory, numerical efficiency and the preservation 
of self-duality. The good news is that the Hermitian matrices will commute when the unitary 
matrices commute, and the mapping will be continuous, even Lipschitz if we are careful. This 
means that if (Hi, H 2 , H 3 ) is bounded away from commuting triples, then (Ui, U 2 ) is bounded 
away from commuting pairs. 

The issue with numerical efficiency is that if the matrices U almost commute, then the 
matrices H also almost commute. However, the norm of the commutator may increase under 
the mapping. By choosing suitable, sufficiently well-behaved, maps, we can avoid a large 
increase in the commutator while still obtaining an efficient numerical procedure. 

The map we use for this purpose is, in terms of periodic coordinates on T 2 and real 
coordinates on S 2 CR 3 , 

7(0,0) = (xi,x 2 ,x 3 ) 

where 

xi = f(<t>) 

%2 = g(4>) + h(4>) cos(2iri9) 
X3 = h(4>) sin(27ri6') 

and /, g and h are continuous, periodic scalar functions that satisfy 

f + 9 2 + h 2 = l 
gh = 0. 

To ensure the mapping T 2 — > S 2 is one-to-one one all but a one-dimensional set, ensuring 
later the correct K-theoiy, we need to be sure / moves just once from —1 to 1 on the set 
where h ^ 0. A choice that works well is 



(5.3) f((f)) = if sin(27T0) + H sin(67T0) + ^ sin(lO7T0) 

9(<P) = 



i-OW -\<o<\ 



i-OW \<o<\ 
-\<e<\ 



h(<f>)-- 

although the earlier mathematical work on the Bott index uses / that was piecewise linear. 
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Figure 5.1. Two sets of functions for mapping the torus the sphere, (a) 
Here we take / to be piecewise linear, (b) Here fi(x) is as in (15. 3p . 




Figure 5.2. Panel d) indicates using /, g and h to map the torus (panel a)) 
close to the sphere. Panels b) and c) show the result of truncating the three 
functions to degrees one and three. 

We can interpret this at the level of almost commuting unitary matrices U\ and Ui in 
several ways. For example 

h(<j>) cos(27ri0) = i (h((f))e- 27rie + h((j))e 27Tie ) , 

but also 

h(4>) cos(2ttz0) = ^ (h((f))e^ ie + e- 2me h(<J)) + h(<J))e 2wie + e 2ni9 h((f))) . 
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As we have many matrix symmetries to worry about, we prefer this second formula, and 
define 



(5.4) 



#i = f(U 2 ) 

H 2 = g(U 2 ) + - A [h(U 2 ), U\) + 1 - {h(U 2 ), U x } 



Computing f(U 2 ), g{U 2 ) and h(U 2 ) is simple enough in theory, as we can diagonalize U 2 
by a unitary W, so 



- g27rixi 



u 2 = w 



o2 nix n 



f(U 2 



w 



f{Xr 



We shall refer mapping the torus to the sphere as in ( 15.41) as mapping by the polynomial 
map as we implement it by calculating nearby polynomials in U\,U\,U 2 ,\j\. There is an 
alternate method that is slower to compute but appears to show phase transitions more 
clearly. We shall refer to it as the logarithmic map and is defines 



(5.5) 



where 



H x = £(U 2 ) 

H 2 = \{m(U 2 )M} + \{m(U 2 ),U 1 } 
H^ % -[m{U 2 )M}-\{rn{U 2 ),U x } 



is discontinuous and m is the continuous function m = \/l — £ 2 . 



Lemma 5.1. With Hi, H 2 and H 3 defined as in (5^4), if the unitaries U\and U 2 satisfy 
II [Ui,U 2 ] || < 5 then 

3 



r=l 



<0{5) 



and 



\[H r ,H s }\\<0(5). 



It is also easy to prove a stability result, to the effect that if we also have V\ and V 2 and 
define the associated three Hermitians Ki, K 2 and K3, then 



where 



\H r -K r \\ < 0(rj) 



r ? = max(||?7 1 -V r 1 ||,||f/ 2 -V 2 | 



The upshot is that we can now use the Bott index of the sphere to determine an index on 
the torus 
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Definition 5.2. If ?7i, are unitary matrices, define H r by (15. 4p and (15. 3p . 

bott(Z7i,E7 2 ) =bott (H X ,H 2 ,H 3 ) 
If C/i, £7 2 are self-dual unitary matrices 

bott {Ui, U 2 ) = bott (Hi, H 2 , H 3 ) . 
Theorem 5.3. Suppose U\ and U 2 are unitary matrices and 

\\[UiM\\<8. 

For small S, the torus index of unitaries U\ and U 2 equals 

Bott (Hi,H 2 ,H 3 ) 

where H r are as above, with either choice for f, g and h, or where the H r are defined by the 
log-method. 

Proof. This is essentially the main theorem in |T7] . The only new claim is that the smoother 
choice of /, g and h can be used. There is a continuous deformation between the two choices 
of functions, keeping the gh = and f 2 + g 2 + h 2 = l relations at all times, so for small 
enough 8 the resulting path of Hermitian matrices will keep a gap in its spectrum at 0. 
Therefore the eigenvalue counts are not able to vary. □ 

Conjecture 5.4. For small 8 the index bott (Hi, H 2 , H 3 ) is the same when /, g, h are defined 
by either the polynomial map of the logarithmic map. 

5.2. Other geometries. There are other ways to describe the torus than by two unitaries. 
For example, we could consider the torus as embedded in 3-dimensional space. We would then 
have matrices X%, X 2 , X 3 obeying appropriate algebraic relations describing this surface. In 
general, we can consider other spaces in this way. For example, we could describe a system 
on a two-dimensional manifold with many handles by imagining this manifold embedded 
in 3 dimensional space. We would then construct three matrices Xi,X 2 ,X 3 again with 
appropriate algebraic relations. We then project these matrices as before. This gives us a 
soft representation of the original manifold. 

For each such space, we can construct a group, called the reduced K , describing possible 
topological obstructions in the GUE case. For example, for a three dimensional torus, we 
find that this group is equal to Z + Z + Z. These three integer invariants are in fact lower 
dimensional invariants, similar to the idea of weak topological insulators studied in the 
translation invariant case [19]. These can be understood as follows. We have three matrices, 
Ui,U 2 ,U 3 , which are approximately unitary and which approximately commute with each 
other. Any pair of them, such as U\ and U 2 can have a nontrivial invariant as in the case of 
the two torus. The physical interpretation is that we simply ignore one of the three directions 
of the torus. We have a Hamiltonian which is local on the three torus, and then we simply 
map the lattice sites to sites on the two torus by ignoring one of thre three coordinates and 
we then construct a Hamiltonian which is local on the two torus. For other manifolds, just 
as in the case of the three-dimensional torus, we may see topological obstructions arising 
from lower dimensions; these will appear in the reduced Kq. 

One way of obtaining just the highest dimensional obstruction, is to map the manifold 
onto S +1 . Thus, we map T 2 onto S 2 , or T 3 onto S 3 , and so on. This generalizes the torus 
to sphere mapping described above. In the case of a GUE system, this leads to nontrivial 
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obstructions only for d even. We will also use this mapping procedure in the case of three- 
dimensional time-reversal invariant insulators described below. We study a system whose 
lattice sites live on a three-dimensional torus. However, after constructing the projected 
position matrices on the torus, we then map these matrices to matrices on a sphere, and 
then study invariants on the sphere. One reason for this is that we do not have a "native" 
torus formula for the index for any system other than the complex case on the two-torus. 
In all other symmetry classes, we only know how to compute the index by mapping to the 
sphere. 

In the special case of the three torus, we can define the needed four almost commuting 
self-dual Hermitians by one of two methods. Suppose Ui, U 2 and U 3 are self-dual almost 
commuting matrices. We first define 

M = f(U 3 ) + ig (U 3 ) + \h (U 3 ) f (U 2 ) + if (U 2 ) h {U 3 ) 
N=h (U 3 ) g (U 2 ) + \h (U 3 ) h (U 2 ) U x + \U x h (U 2 ) h (U 3 ) 

which gives self-dual matrices that are almost normal and almost commute. From there we 
set 

H t = \ (M f + M) 
YL 2 = \ (M f - M) 
H 3 = \ (JV+ + N) 
# 4 = f (iVt - N) 

This is our polynomial map. 

For the logarithmic map, we define 

M = £(U 3 ) + \m (Us) f (U 2 ) + |/ (U 2 ) h (U 3 ) 

N=h (U 3 ) g (U 2 ) + \h (U 3 ) h (U 2 ) U x + \U x h (U 2 ) h (U 3 ) 

and then define the H r as above. 

5.3. Relation of Matrix Invariant to Hall Conductance. We now show that the matrix 
invariant of the matrices U\ , U 2 that we compute is the same as the usual Chern number 
invariant. This will in fact provide a simple proof of Hall conductance quantization for free 
fermion systems. 
Consider 

(5.6) ta(log(tfiE/ 2 E#E4)). 

We have shown that this is equal to r+2mm, for some integer m, with r = log( | (iet(f/ 1 ) | 2 |c?et(f/ 2 ) | 2 )- 

Lemma 5.5. Assume U±,U 2 are obtained from a free fermion system on a torus topology. 
Let the lattice be a square lattice, of size L-by-L, with L 2 lattice sites. Let the free fermion 
Hamiltonian have hopping distance bounded above (uniformly in L) by a constant, and spec- 
tral gap bounded below (uniformly in L) by a constant. Then, \\[P, Ui}\\, \\[P,U 2 }\\ < 0(1/L), 
and || [Ux,U 2 \ \\ < 0(l/L 2 ) and \\UxU\ - \\U 2 U\ - I\\ < 0(1/L 2 ). 



Proof. This is a minor variation of lemma 5.1 in [TJ. 



□ 



U x XJ 2 Xj\Xj\-I <0(1/L 2 ), 
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This implies that 

(5.7) 
so 

(5.8) || log^iWM) - (UxU 2 U\u\ - ij || < C(l/L 4 ). 
Thus, 

(5.9) Im(tr(lW2t/M - m)) < 0(1/L 2 ), 

using the fact that the dimension of the Hilbert space is L 2 and the trace of an operator is 
bounded by the dimension of the space times the operator norm of that operator. 
Note that 

(5.10) rm(tr([/if/ 2 f/M)) 

= Im(tr([fA, U 2 \U\UI)) + Imitr^UiUM)) 

= lm(tT([U u U 2 ]uM)), 
where we used the fact that U a U\ is Hermitian, and given any two Hermitian matrices A, B, 
we have Im(tr(AB)) = 0. 

For notational convenience, let us define 



(5.11) 

(5.12) 
so that 

(5.13) 



(5.14) 

Then, 
(5.15) 



U = Pexp(iQ 1 )P, 
V = Pexp(iQ 2 )P, 



U 



V 





U x 



u 2 



lm(ti([U u U 2 }UlUl)) = Im(tr([C/,y][/ t y t )). 



Define the current operators J\ , J 2 by 



J a = i 



.exp(i6 a )if exp(-iO a ) - exp(-i6 a )if exp(i6 a ) 



(5.16) 
Then 

(5.17) J a = [exp(*0 a ), H] exp(-20 a ) + 0(1/ L 2 ). 

The first equation for the current operators is Hermitian; the second is slightly more conve- 
nient for later use. Let ^(1 — P)J a P S j(t) = exp(iHt)(l — P)J a Pexp(—iHt), following the 
usual Heisenberg evolution of operators. Then, 

(5.18) im + y dtexp(et)(l - P)J a (t)P = lim J dt exp(rf) ((1 - P) J a P) (t) 

= _j(i _ p) exp (^e a )Pexp(-20 a ) + 0(1/L 2 ) 
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To obtain the above equation, we first use use Eq. (I5.17P to approximate J a . We then 
set (1 - P)[exp(tO a ),H]exp(-tQ a )P = (1 - P)[exp(iQ a ),H](l - P)exp(-tO a )P + (1 - 
P)[exp(i6 a ), P]Pexp(— iO a )P. The first term is bounded by 0{1/L 2 ) since we can bound 
the norms of [exp(i0 a ), H] and (1 — P) exp(— iQ a )P both by 0(1/ ' L). We then approximate 
(l-P)[exp(te a ),H]Pexp(-tQ a )P=[(l-P)exp(tQ a )Pexp(-tQ a )P,H]+0(l/L 2 ). Finally 
we use that lim e _^ + /_ dtexp(et)[0, H](t) = —iO for any operator O which has vanishing 
matrix elements between degenerate eigenstates of H . Similarly 

r° 

(5.19) lim / dtexp{et)PJ a (t)(l- P) = -iPexp(i9 a )(l - P) exp(-i0 a ) + 0{1/L 2 ) 

-2expHe a )Pexp(20 a )(l - P) + 0(1/L 2 ). 
Let 

(5.20) 0+ = (1 - P) exp(ze 1 )Pexp(-26 a ), 
and 

(5.21) 0~ = exp(-ze a )Pexp(26 a )(l - P). 

Using the fact that [exp(z©i), exp(z0 2 )] = 0, we have P[exp(i6i), exp(i6 2 )]P = 0. How- 
ever P[exp(«ei),exp(ie 2 )]P = [U,V] + Pexp(i0i)(l - P) exp(i0 2 )P - Pexp(i6 2 )(l - 
P) exp(i9i)P. So, [U, V] = Pexp(z6 2 )(l - P) exp(z9i)P - Pexp(i9i)(l - P) exp(i6 2 )P. 
Thus, 

(5.22) lm{tr{[U,V}U j V j )). 

= Im(tr(Pexp(i6 2 )(l -P)exp(i6i)Pexp(-iei)Pexp(-ie 2 ))) 
-Im(tr(Pexp(i6i)(l - P) exp(z0 2 )Pexp(-i6i)Pexp(-ie 2 ))). 

Since ||[P, exp(i6i)]|| and || [P, exp(i0 2 )]|| are both bounded by 0(1/L), the term inside the 
first trace in the above equation is equal to Pexp(z0 2 )(l— P) exp(z9i)Pexp(— zGi) exp(— z0 2 )P, 
up to an error in operator norm which is of order 1/L 3 . Since the dimension of the Hilbert 
space is equal to L 2 , this means that the first trace is equal to Im(tr(Pexp(i0 2 )(l — 
P) exp(iei)PexpH9i) exp(-z9 2 ))) + 0(l/L) = -Im(tr (OfO^) + 0(1/L). Applying the 
same argument to the second trace in the above equation, 

(5.23) lm{ii{\U,V}U ] V ] )). 

= Im(tr(0+0 2 -0+0 2 )) + 0(l/L). 
Using the result (15.181) this becomes 

/o rO 
dtiexp(eiti) lim / dt 2 exp(eit 2 )Im(tr(J 2 (t 2 )(l - P) Ji(ti)P - Ji(ti)(l - P) J 2 (t 2 )P). 
.00 ^2^0+ J_ 00 

Recognizing this as the Kubo formula for the Hall conductance, we find two results. 

Lemma 5.6. On a torus, the matrix invariant is equal to the Hall conductance, up to 
0(1/L). 

This has the corollary: 

Corollary 5.7. Consider a free fermion system on a torus topology. Let the lattice be a 
square lattice, of size L-by-L, with L 2 lattice sites. Let the free fermion Hamiltonian have 
hopping distance bounded above (uniformly in L) by a constant, and spectral gap bounded 
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below (uniformly in L) by a constant. Then, the Hall conductance is within 0(1/ L) of an 
integer. 

While better bounds are known and are valid for interacting systems [30]. this seems a 
simple way to show quantization of Hall conductance for non-interacting systems. See also 
[32] for the non-commutative geometry approach to this problem and |33j for a review of 
non-commutative geometry and topological insulators. 

This result also implies that 

Lemma 5.8. Assume U\,U2 are obtained from a free fermion system on a torus topology. 
Let the lattice be a square lattice, of size L-by-L, with L 2 lattice sites. Let the free fermion 
Hamiltonian have hopping distance bounded above ( uniformly in L) by a constant, and spec- 
tral gap bounded below (uniformly in L) by a constant. Then, for all sufficiently large L, the 
matrix invariant agrees with the free fermion invariant in both the GUE and GSE universality 
classes. 

Proof. We have just shown that the Chern number invariant, which is the same as the free 
fermion invariant, agrees with the matrix invariant in the GUE case[8]. Consider the GSE 
case. We use a similar argument as in the paragraph near Eq. (13. ip . adapted to the Z 2 case. 

Our first step is to construct a free fermion Hamiltonian in the GSE class with free fermion 
invariant and matrix invariant both equal to —1. Let Hque be a free fermion Hamiltonian 
in the GUE class, with Chern number +1. Define a free fermi Hamiltonian Hgse by taking 
two copies of Hque- 



where the overline denotes the complex conjugation. By identifying the two copies with spin 
up and down, the Hamiltonian Hgse is indeed time- reversal invariant and has free fermion 
invariant —1. Since Hque has odd Chern number, the corresponding band projected position 
matrices of ifc-i/shave odd matrix invariant. The band projected position matrices of Hgse 
are obtained by doubling the band projected matrices of Hque] that is, given matrices H r 
from Hamiltonian Hque, consider the matrices H' r defined by 



Given that H r have odd matrix invariant, the doubled matrices have Z 2 invariant equal to 
— 1, as desired (see theorem (12. 8p . 

Consider any free fermi Hamiltonian H in the GSE universality class. If its free fermion 
invariant is equal to +1, then H is equivalent, up to addition of trivial degrees of freedom, 
to a Hamiltonian H' with localized Wannier functions (here we use a result of Kitaev in[8], 
though the details of that proof are not yet published). This implies that, for sufficiently 
large L, the matrix invariant of H' is also equal to +1. This implies that the matrix invariant 
of H plus trivial degrees of freedom is equal to +1, which implies that the matrix invariant 
of H is equal to +1. Therefore, if the free fermion invariant is equal to +1, then the matrix 
invariant is also equal to +1 for sufficiently large L. 

Suppose instead H has free fermion invariant equal to —1. Consider the Hamiltonian 
H (B Hgse- This has free fermion invariant equal to +1. Following the same argument as 
above, H © Hgse has matrix invariant equal to +1. Using the fact that the matrix invariant 



(5.24) 




(5.25) 
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of if © Hqse is the product (using the Z 2 group multiplication rule) of the matrix invariant 
of H with that of Hqsei the matrix invariant of H is equal to —1. □ 

6. Chiral Classes 

The discussion above is entirely concerned with the classical universality classes. In this 
section, we consider the three chiral symmetry classes (unitary, orthogonal, and symplectic). 

We begin with the complex case. In the complex case, the chiral class consists of Hamilto- 
nians on a bipartite lattice, with non-zero hopping only between sites on different sublattices. 
Such Hamiltonians can be written as 

«"> (J. o) ■ 

where the two blocks correspond to the odd and even sublattices, respectively. In the complex 
case, chiral Hamiltonians have topological obstructions in odd dimensions, while non-chiral 
Hamiltonians have topological obstructions in even dimensions [8]. 

However, the topological obstructions for chiral Hamiltonian do not arise from obstructions 
to construct localized Wannier functions, unlike the other problems we study. Consider a 
simple example in one-dimension. We have a system with N sites arranged on a ring, with 
N even. Let Hamiltonian ifo be 

(6.2) i2o = +/i.c, 

i=2k 

and let Hi be 

(6.3) H i= Yl ^i+i + h.c. 

i=2k+l 

These two Hamiltonians are both gapped, but they are in different topological phases: there 
is no way to find a continuous path connecting these Hamiltonians while maintaining locality 
and while keeping the gap larger than 0(1/ N). However, both of these Hamiltonians have 
localized Wannier functions. For the Hamiltonian H , the localized Wannier functions for 
the occupied states (we fix Ep = for all chiral Hamiltonians) are given by N/2 vectors, v , 
where v k has entries (v k )i = l/y/2 if i = 2k, (v k ) i+1 = —l/y/2 if i = 2k + 1, and (v k )i = 
otherwise. That is, these vectors are localized on sites i,i + l for even i. For Hamiltonian 
Hi, the Wannier functions are localized on sites i,i + l for odd i. Thus, the difference in the 
phases does not have to do with one Hamiltonian having localized Wannier functions and 
the other Hamiltonian not having localized Wannier functions. 

However, we can still quantify topological obstructions in chiral systems using almost 
commuting matrices in a different way. We consider the case without time reversal symmetry 
first. Let if be a chiral Hamiltonian. Spectrally flatten H to define a Hamiltonian H' whose 
eigenvalues are all equal to ±1. Using P as the projector onto negative eigenvalues of H, we 
have the relation 

(6.4) H' — 1 — 2P. 

The Hamiltonian H' is still chiral. Suppose H has a gap in the spectrum near so that 
all eigenvalues of H are greater than AE in absolute value; then H' is still local so that if 
entries of H decay superpolynomially in the distance between sites then so do entries of H' 
and if entries of H decay exponentially then so do entries of H' (this can be shown using 
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the same techniques used to prove exponential decay of correlation functions [31]). Further, 
since H' is an odd function of H, then H' is chiral given that H is chiral. 

For definiteness, let us work on a d- dimensional torus. We define position matrices exp(i@j) 
for % = 1, d for sites on the even lattice exactly as we did in the non-chiral cases. We then 
pair sites in the lattice, choosing pairs of neighboring two sites on opposite sublattices and 
considering them to be a "pair" of sites (we assume the total number of sites in the lattice 
is even; if it is odd, then there are zero modes in H). For example, in the one- dimensional 
system above, we may choose to pair sites 1 and 2, sites 3 and 4, and so on. Since H' is 
chiral, we can write H' in the form f)6.ip . We use this pairing in writing H' in the form (16. 
if there are a total of N sites in the lattice, we choose the i + iV/2-th basis vector to be the 
pair of the i-th basis vector. The first N/2 basis vectors are in one sublattice and the last 
N/2 are in the other. Since H' 2 = I, the matrix A is a unitary matrix. 

Thus, we have constructed d + 1 different unitary matrices: d of these are the matrices 
exp(iOj), while the d + 1-st is the matrix A. These matrices almost commute with each 
other, since [exp(iGj), exp(zOj)] = and using locality properties of P we can bound the 
operator norm of the commutator [exp(iOj), A]. Non-trivial topological obstructions can 
exist for d + 1 unitaries. For example, consider the Hamiltonian Hq, Hi given above. This is 
a system on a 1-torus, so we have 2 unitaries, both of size N/2-by-N/2. One unitary is the 
diagonal matrix 



(\ 



(6.5) exp(i6) = 

The other unitary is equal to 
(6.6) 

for H , but it is equal to 
(6.7) 



exp(i27r/(iV/2)) 



/l 



A, 



exp(i4vr/(iV/2)) 



\ 



7 



An 



'0 1 
1 



...1 



for Hi. The matrices exp(zO) and A exactly commute. The matrices exp(z0) and Ai 
almost commute (the operator norm of the commutator is of order 1/iV), but cannot be 
approximated by exactly commuting matrices, as can be seen by computing the invariant 
m from Eq. (15.21) . In fact, these two matrices exp(iO) and A\ are a previously considered 
example of a pair of almost commuting unitaries which cannot be approximated by exactly 
commuting unitaries [18]. 

Two almost commuting unitaries are characterized by an integer invariant as described in 
the previous section on the torus, corresponding to the integer invariant known to describe 
chiral systems in one dimension. Note that this invariant of chiral systems provides an ob- 
struction to connecting two Hamiltonians by a path of gapped, local Hamiltonians. Suppose 
H , Hi are gapped, local Hamiltonians, connected by a smooth path H s of gapped local 
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Hamiltonians. Since H s is gapped and local, the corresponding spectrally flattened Hamil- 
tonian is local, and so the unitary matrix A s is local. Thus, the matrix A s approximately 
commutes with exp(zB) for all s. If the commutator [A s , exp(zO)] is sufficiently small, then 
the integer invariant described above is indeed invariant under small changes in A s . Thus, 
if the gap is sufficiently large, the invariant does not change along the path of Hamiltonians 
and so is the same for H Q and Hi. 

The system of d + 1 unitaries describes a soft d + 1-torus. In the torus case, we have 
obstructions for all d > 1, due to the possibility of lower dimensional obstructions, just as in 
the case of weak topological insulators discussed in the non-chiral case, however the highest 
dimensional obstruction occurs only for d odd. We can repeat the exercise on a sphere, 
instead of a torus. In this case we obtain d + 1 Hermitians Hi obeying the requirement 
that Yli H? = I-i an d one unitary A which almost commutes with the H^. This gives a soft 
S d x S 1 . The reduced Kq of S d x S 1 is Z in all dimensions. For d even, our system is always 
in the trivial case, because the topological obstructions for this space in even d occur only 
if the d + 1 matrices describing S d have a topological obstruction, and in our case these 
d + 1 matrices commute exactly. So, we see integer obstructions in odd dimensions and no 
obstructions in even dimensions. 

One can also consider the chiral real and self-dual classes in this manner. In this case, the 
matrices A are orthogonal or symplectic matrices, respectively. 

The same mathematical problems are present in the case of sublattice symmetry as in the 
non-chiral case. We need to show that the index obstructions we have obtained are the only 
obstructions, and to construct explicit examples of sequences of almost commuting unitaries 
displaying the different obstructions. Finally, we need to identify this invariant with other 
known invariants. 

7. Numerical Results in Three Dimensions 

We now describe our numerical results on a three dimensional time reversal invariant 
topological insulator. We considered the system on a three dimensional torus, using the 
polynomial map described previously to map to the sphere. 

7.1. Three Dimensional Hamiltonian. In previous numerical work in two dimensions [9], 
we used the model of [31] which includes coupling between up and down spin components 
due to breaking of bulk inversion symmetry, plus an additional coupling to disorder. The 
model we study in three dimensions consists of a Hamiltonian 

(7.1) U = U + V, 

where V represents on-site disorder described below. The Hamiltonian "H representing the 
system without disorder is defined on a three dimensional torus. Each site has four different 
states, corresponding to two different bands and to two different spins. We introduce Pauli 
spin matrices cr b x a ™ d to describe band desgrees of freedom and cr x P y n z to describe spin degrees 
of freedom. We set 

(7.2) H = A{id x a b x and af n + id y a b y and + id z a band a s x pin ) + (M + Bd 2 )a band , 

where A, B, M are numerical parameters, and d Xiy ^ z is short-hand notation for a lattice 
derivative. That is, d x denotes a matrix whose matrix elements between sites j and k is 
equal to 1 if site j is one lattice site away from site k in the x-direction, equal to —1 if 
site j is one lattice site away from site k in the — x-direction, and otherwise. Similarly, 
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Figure 7.1. Plot of variance for L = 4, 6, 8, 10, 12. 



<9 2 is a lattice second partial derivative: its matrix element between sites j and k is equal 
to —6 if j = k, equal to +1 if j and k are nearest neighbors, and otherwise. We chose 
A = 1, B = —1, M = —2. These parameters were chosen to obtain a topologically nontrivial 
phase in the absence of disorder, with the relation between B and M chosen to cancel the 
leading irrelevant term k 2 in a continuum treatment of the problem. This Hamiltonian 
appears in [25] . 

The disorder term V is a diagonal matrix, independent of spin and band index. The 
disorder on a given site was chosen randomly from the interval [—4,4]. This is sufficiently 
strong to close the gap in the middle of the spectrum. The spectrum of H is symmmetric 
about zero energy, and the distribution from which we draw V is also symmetric about zero 
energy. The statistical properties of the spectrum of H are symmetric about zero energy. 

We checked the localization properties of the eigenvalues as follows. For each disorder 
realization, for each eigenfunction, if), we computed the variance in sin(^), cos(^) for each 
of the three angles 6*i, 9 2 , O3 on the three torus. We averaged over angles, computing 

1 3 

( 7 - 3 ) 3 (Z)^' M^V) + cos(0i) 2 V> - (V, sin(^> 2 - (V>, cos( W> 2 ) 

i=l 
1 3 

= 1-3 (I>> sin ( W> 2 + cos(^> 2 ) . 

1=1 

. We plot this variance as a function of energy in Fig. (17. ip . One can see crossings in the 
variance as a function of system size at E « ±0.65 and E « ±10.6. For —10.6 < E < —0.65 
and 0.65 < E < 10.6, the variance increases as a function of system size, indicating that 
the eigenfunctions are delocalized in this energy range, while outside this range the variance 
decreases, indicating that the eigenfunctions are localized. Let E[ oc ps —10.6 and E[ oc w —0.6 
denote the localizations of these two localization transitions. 

When we turn to the index as a function of energy, we encounter a surprise. The expected 
behavior (which was seen in two dimensional systems in [TJ), is that the delocalized region of 
energies should precisely coincide with the region in which the index fluctuates from sample 
to sample, while in the localized region the index should be the same in all samples. The 
surprise we encounter in three dimensions is that there is a range of energies over which the 
index is equal to ±1 in all samples, but the eigenfunctions are still delocalized. In Fig. (17. ip 
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we plot the index computed using a polynomial mapping from torus to sphere. For Ep less 
than approximately —4, the index is equal to +1 for almost all samples, and it seems likely 
that in the thermodynamic limit, the index will be equal to +1 for all samples for Ep < E c 
for some critical energy E c which is roughly —5. Conversely, it also seems likely that in the 
thermodynamic limit the index will be equal to —1 for all samples for E' c < Ep < 0, for some 
critical energy E' c . This energy E' c appears to be at around —0.6. 

Since for some small fraction of samples the index was not equal to +1 for Ep = —4.9 
(the most negative Ep plotted in the figure), we ran additional checks on a larger number of 
samples with Ep = —5, 6, 7, 8, 9, 10. We are able to compute the index more rapidly for the 
samples with more negative Ep since the number of occupied states decreases. We studied 
« 2700 samples for L = 10 and « 1800 samples for L = 12. We found that for L = 10, the 
average index was equal to 0.998... for Ep = —7 and was equal to 1 for Ep < —8, and for 
L = 12 the average index was equal to 0.992... for Ep = —6 and was equal to 1 for E F < —7. 
This provides strong numerical evidence that in the thermodynamic limit the index is indeed 
equal to 1 for Ep < E c « —5. Perhaps E c is slightly smaller than —5, but certainly the 
numerical evidence indicates that it does not coincide with the localization transition. 

Thus, the critical energy E' c where the index becomes equal to —1 appears to coincide 
with the transition to localized states at E' loc , while the critical energy E c appears to be 
in the middle of the delocalized region. This is surprising, but does not contradict any of 
the known properties of the index. Note that the index must be trivial (+1) if all of the 
occupied states are localized, because in this event we know that the system has localized 
Wannier functions: the energy eigenfunctions themselves supply a basis of localized Wannier 
functions. Hence, the index must be equal to —1 for Ep < Ei oc . 

Similarly, the fact that the index is fluctuating from sample-to-sample for E c < Ep < E' c 
implies that the system must be in a delocalized phase in this region. To state this claim 
in a more mathematically precise fashion, the fact that the index fluctuates from sample-to- 
sample and monotonically decreases with Ep, implies that for any small interval of energies, 
[E, E + dE], with E c < E < E' c , there is some non-zero probability that a given sample will 
have a transition in index from +1 to —1 in that interval, which means that there is some 
non-zero probability that a state in that interval will be delocalized. 

The surprising thing, then, is that the index can be equal to +1 for all samples even 
when the system is delocalized. However, perhaps this should also not be surprising. Recall 
the phenomenon in the two-dimensional casejl]. There, the index always fluctuated from 
sample-to-sample in a delocalized phase. However, the average index varied smoothly as a 
function of energy in the delocalized phase. The average was equal to +1 exactly at the 
lower critical energy E c , and the average was close to +1 cose to the the lower critical energy 
E c . Thus, the system was delocalized and yet the index was be equal to +1 in almost all 
samples. Thus, perhaps it should not be surprising that one can have a system in which the 
index is equal to +1 in all samples. 

The question that this raises is how to interpret the index in a delocalized phase. In 
a localized phase, the index quantifies topological properties of the system, and we have 
proven that it does not change under small change in the Hamiltonian (if the Hamiltonian 
is localized, small changes in the Hamiltonian only lead to small changes in the projector 
P and hence to only small changes in H r and in B(Hi, Hd+i)). In the delocalized phase, 
we do not have a good interpretation of our index. However, it has been pointed out to 
us|24j that the index is well-defined so long as the entries of P drop off as a sufficiently rapid 
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Figure 7.2. Plot of average index for L = 6,8,10,12 (cir- 
cle, square, diamond, star), using polynomial map from the three-torus to 
the three-sphere. Each data point is an average of 1700,1400,600,400 
samples, respectiv 



G-Q L=6 
B-B L=8 

L=12 




Figure 7.3. Plot of average index for L = 6,8,10,12 (cir- 
cle, square, diamond, star), using logarithmic map from torus to sphere. 
Each data point is an average of 1300, 1000, 300, 270 samples, respectively. 



power of distance (the particular power depends on spatial dimension) and hence one can 
have a delocalized phase with a well-defined index. However, the calculation of the index 
reveals an unsuspected critical point in the system in three dimensions, a transition at E c 
which appears not to have any signature in the localization properties of the system. Thus, 
it would be interesting to understand this transition from a field theoretic point of view. 

7.2. Particle-Hole Symmetry and Index. In the calculation in this section, we computed 
the index only for E F < 0. The reason to do this is that the dimension of the matrix B is 
equal to twice the number of occupied states. For E F < 0, this number of occupied states 
is less than half the dimension of the Hamiltonian, so that the dimension of B is at most 
equal to the dimension of the Hamiltonian. This restriction on the size of the matrix B is 
important in making it practically possible to consider large system sizes. We implemented 
the algorithm by first diagonalizing the Hamiltonian. Then, we computed the index for a 
sequence of values of E F , starting at E F = and reducing by 0.1 at each step. As E F gets 
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Figure 7.4. Plot of average index for L = 10, comparing logarithmic and 
polynomial maps (circle, square). 

more negative at every step, the dimension of B decreases, so that the calculation of the 
index is faster for more negative values of Ep. 

However, one may wonder what happens for Ep > 0. Is there new physics? In this 
subsection, we discuss the implications of particle-hole symmetry for Ep > and argue that 
we can in fact obtain all the desired information from Ep < 0. Assume that we have some 
system with a given Ep. Let P be the projector onto the space of occupied states. Then, 
1 — P is the projector onto the space of unoccupied states. We can form the matrices 



where (1 — P) is the projector onto the space of empty states. We need the following lemma: 

Lemma 7.1. For any symmetry class from GUE,GOE,GSE and any dimension, there is 
a numeric constant c > such that the following holds. Assume that the commutator 
\\[X a ,P}\\ < c for all a and that Yla-^a = I an ^ [X a ,X b ] = 0. Compute the the index 
from the matrices H' r and call this I'. Compute the index from the matrices H r and call this 
I. Then, the product of I with V , using the appropriate group multiplication law (either Z or 
1*2), is equal to the trivial index. Equivalent, I = I' in the Z2 case, while I = n and V = —n 
for some integer n in the Z case. 

Proof. We need two properties of the index, proven elsewherepQ. First, we use the fact that 
under direct sum of matrices, the index combines using multiplication as discussed above 
Eq. 13.11 in the Z case. That is, let I({H r }) denote the index computed from a given set of 
matrices H r . Given any matrices H r and H' r , we have 



So, to show the claim it suffices to show that I({H r © H' r }) is trivial. However, we have 





(7.7) 




(7.8) 



H T ®H' r = X a - PX a (l - P) - (1 - P)X a P. 
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If ||[X a ,P]|| < c, then \\PX a {l - P)\\ < a So, \\H r © H' r - X a \\ < 2c. The matrices X a 
have trivial index. The second property of the index we need is that one can prove a lower 
bound on the distance in operator norm from the X a to the nearest set of matrices with 
nontrivial index (proven in [TJ for the case of d = 2 in two different symmetry classes, and 
stated without proof here for other cases since the proof is essentially the same for all cases). 
Thus, if c is smaller than this lower bound, the desired result follows. □ 

Then, suppose the commutator ||[X a ,P]|| is indeed sufficiently small, as would hold in 
the localized regime. Then, using the lemma, we can either compute the index by consider 
the space of occupied states or the space of empty states, choosing whichever space has 
the smaller dimension to make the numerics simpler. In the case of the three dimensional 
system given here, for Ep > 0, we would choose to compute the index from the space of 
empty states. However, using symmetry properties of the Hamiltonian and the disorder 
distribution, one can show that the statistical properties of the index computed from the 
empty states at energy Ep are identical to the statistical properties of the index computed 
from the occupied states at energy — Ep. Thus, assuming the commutator is sufficiently 
small, indeed we do not need to consider the properties of the system with Ep > 0, at 
least in the localized regime. However, on physical grounds we believe that even outside the 
localized regime the indices I({H r }) and I({H' r }) are related in the same way, so long as the 
projector P is obtained from a local Hamitonian. 

8. K-Theory 

We have described many of our invariants in terms of linear algebra and matrix functions 
(functional calculus to pure mathematicians). This makes finding efficient algorithms, but 
hides where these invariants are coming from. 

In finite models, the position observables are bounded Hermitian operators that commute. 
The spectral theorem applies. This immediately brings into the discussion a commutative 
algebra of functions on a space. Let us adopt the notation 

C(X) = { f : X ->■ C | / is continuous } 

whenever X is a compact Hausdorff space. Given commuting n-by-n Hermitian matrices 
Hi, ... , Hk, of norm at most one, the spectral theorem states that there is a unitary W so 
that H r = WD r W^ with D r diagonal, 

D r = diag(A r ,i, . . . , A r ,„). 

The requirement that the H r be contractions (norm at most one) translates to the restriction 
— 1 < \ r j < 1. For / in C([— 1, l] k ) we now are able to define 

/(#!, ...,H k ) = Wdiag(A r ,i, . . . , X r ,i)Wl 

On the one hand, this is a nice form of functional calculus (multivariate matrix functions). 
One the other, we get a homomorphism of algebras 

(8.1) $:C([-l,lf)^M fc (C) 

defined by 

*(/) = f(Hi, ■ • .,-Hfc), 
where A C M n (C) is the unital complex algebra generated by the H r . 
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If the H r satisfy some equation, that equation allows us to cut down the space [— 1, l] fc to 
something more interesting. For example, if 

H I + Hi = Hi + Hl = I 

then we conclude 

\i + ^l,e — ^l,e + ^l,e — 1 

and so we get a homomorphism 

* : C(T 2 ) -> M fc (C). 

We often think of this in terms of two commuting unitaries, U — Hi + iH 2 and V = H 3 + 
iH^. From two commuting unitary matrices U and V in M/-(C) we derive a homomorphism 

: C(T 2 ) M fe (C). 

For Laurent polynomials \l/ operates as expected. For example 

^ (cos(2tt^i) + e^e- 1 ^ 2 ) = \{U + U ] ) + UV~ l . 

The algebras A, M n (C) and C([— 1, l] fe ) are examples of complex C*-algebras. The cate- 
gory of complex C*-algebras gives us enough room to study algebras of matrices generated 
by almost commuting Hermitian matrices. If we wish to keep track of TR-symmetry via 
self-dual matrices, we need to switch to the study of real C*-algebras. These can be thought 
as complex C*-algebras that carry an extra operation, similar to the matrix transpose or 
matrix dual. 

8.1. Complex C*-Algebras. C*-algebras are algebras over the complex scalars that are 
equipped with a norm and an involution a;4i*. Only in the context of real C*-algebra does 
one generally specify complex C*-algebra. 

There are axioms relating the algebraic and norm structures, c.f. [20]; the idea is that 
the norm is modeled on the operator norm, while the involution is modeled on the adjoint 
of a complex matrix or of an operator on Hilbert space. The adjoint is written as A* in 
mathematics literature and as A* in physics literature. We will use * for the involution in an 
abstract C*-algebra, and f when dealing with finite matrices or actual operators on Hilbert 
space (for example, when considering matrices above, we used f for the adjoint throughout). 

The most basic mappings between C*-algebras are the *-homomorphisms. A *-homomorphism 
if : A — > B must be linear and satisfy the axioms 

(f(ab) = (f(a)if(b) 

and 

(p(a*) = (p(a)*. 

In examples the involution should be represented appropriately. For example, the involution 
in C(X) is conjugation, the value of / at x is f(x). The map $ in ( 18. II) is a ^-homomorphism 
because 

Hf) = Hf) j 

Continuity of a *-homomorphisms is automatic, as is the fact that when (p is ^-homomorphism 
we have ||<^(a)|| < ||a||. 
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8.2. Complex K . There are two groups K Q (A) and K\{A) that make up the i^-theory of 
a complex C*-algebra A. A straight-forward, common construction of K (A) uses projectors 
(a.k.a. projections, a.k.a. self-adjoint idempotents) and homotopy. The simplest construc- 
tion of Ki(A) uses homotopy and unitaries. The ability to build up matrices is used to gain 
a commutative group operation. 

Assume A to be a unital C*-Algebra. Were A not unital, we would need to consider 
the so-call reduced K group. We are glossing over such details in this short summary of 
.fT-theory in the complex case. See [21] for details. An alternate picture can be used, where 
invertible Hermitian elements replace projectors. Given x invertible with x* = x we can 
associate the projection 

p=\ + |polar(x). 

When x is self-adjoint, polar(rc) will be self adjoint, which forces p to have spectrum in {0, 1} 
and to be a projector. 
Let 

GL° dd (A) = GL n (A) = {xe M n (A) \ x~l exists} 

and 

GL^(A) = {xe Gh n {A) | x* = x} 

To build Kq(A), one considers the union over all even k of the sets GL^(A) and then forms 
equivalence classes denoted [x\. The equivalence relation is defined by requiring 

x ~ y [x] = [y] 

where ~ is to indicate homotopy within the set of invertible Hermitian elements and 

x 

1 

-1 

The group Kq(A) consists of all the classes [x] and formal differences of classes [x] — [y]. 
The addition is determined by the rule 

w + m =(;;), 

while the inverse is formal, 

-([*] - iv]) = [y] - N- 

An example of a K Q class when A is the commutative C*-algebra A = C(X) is that 
associated to a vector bundle. If that vector bundle is a sub-bundle of the trivial bundle 
C'xl then there is a continuously varying projector p x so that the fiber over x is p x (C k ). 
The class in Kq(C(X)) associated to this vector bundle is [2p — 1] where p is the function 
p(x) = p x which is a projector in 

M k (C(X)) = C(X,M k (C)). 

A more basic example is Kq(C). In this case we get an isomorphism 

(8.2) Sig,:# (M„(C))-^Z 
determined by the signature [X] !->■ Sig (X) as in Definition 12.141 



= [*]. 
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8.3. Complex K 1 . We can build Ki(A) with a little less work. First consider the union over 
all k of the sets GL k (A) and form equivalence classes denoted [x]. The equivalence relation 
is defined by requiring 

x ~ y =>• [x] = [y] 
where ~ is to indicate homotopy within the set of invertible elements and 

= [x]. 

The inverse is 

-[*] = [x-i] 

which works because for x and y invertible elements in the same M k (A) the addition has 
two equivalent formulas 

M + \y] = = \. x v\ 

as is shown here: 

(;;)-(S!)(!ir(!!)(!i)~(7!) 

All invertibles in M fe (C) are homotopic — just diagonalize an invertible matrix and push 
the eigenvalues to 1 — so 

Ki{C) = {0} 

where [1] = 0. This stops looking strange after a while. 

Without proof, we claim Ki(C(T 2 )) = Z 2 with the two generators being represented by 
the unitaries u and v in C7(T 2 ) given by u(e 2nie \ e 2ni6 ' 2 ) = e 27Ti01 and v(e 2 ' Ki9 \ e 2nie2 ) = e 2nid2 . 

8.4. Pushing forward /^-theory classes. If we have a *-homomorphism Lp : A — > B be- 
tween C*-algebras, it extends via componentwise action to *-homomorphisms ip k : M.j,(A) — > 
M k (B). These send GL k (A) to GL k (B) and GL™(A) to GL^(S) and so induces two group 
homomorphisms 

<p, : K d (A) -+ K d (B). 

With the examples at hand, the C(X) and M„(C), this is tremendously boring. It is only 
by considering something less rigid than a *-homomorphism that something interesting can 
occur. This dichotomy is the basis for our using i^-theory as a means to distinguish trivial 
an non-trivial phases of models of topological insulators. 

We now give an exact statement about *-homomorphisms (p : C(C) —¥ M„(C). The proof 
uses the fact that homotopic homomorphisms induce the same map on if-theory. By a 
homotopy of * -homomorphisms we mean (p t : A — > B where each (p t is a *-homomorphism 
and with a fixed in A, the map t i— > tpt{a) is continuous. In the following, since C(X) is 
involved, it is more natural to think of Kq in terms of projectors. For projectors in Mfc(C) 
their class in K is determined by the trace. 

Theorem 8.1. Suppose X is a pathwise connected compact metric space. Define I(X) as 
subset of K (C(X)) all elements of the form [p] — [q] where 

Tr(p(x)) = Tr(q(x)) 

for all x. For any unital *-homomorphism 

V :C(X)^M k (C) 
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the induced map 



sends all of I(X) to zero. 



if, : K {C{X)) -> A (M fc (C)) 



Proof. The representation theory of C(X) is well known. By decomposing (p into irreducible 
representations we find a unitary W and points x\, . . . , x k in X so that 



<p{f) = w 



f(x 2 ) 



We can deform W through unitaries to the identity, and deform the Xj along paths to x\ 
and so, as far at A-theory is concerned, ip might as well be the *-homomorphism if> where 

1>(f) = f{x x )L 

Then <p n (p) = p(x\) <S> I and y? n (g) = q{x\) ® I so 

^(fnip)) = A;Tr(p(xi)) = kTr(q(x 1 )) = Tr((p n (q)). 

We are done, by ([821). D 

If we replace *-homomorphisms by mappings that are "almost multiplicative" but otherwise 
like homomorphisms, we can get richer induced maps on i^-theory. This is the basis for 
.E-theory as introduced by Connes and Higson |28j. We prefer the simpler approach of 
defining C(X) by relations and consider our "soft representations" as replacements for *- 
homomorphisms from C(X) to C*-algebras of matrices. These soft representations can 
induce mappings at the level of i^-theory that actual ^-homomorphisms cannot. 

The induced mappings work roughly as follows. A detailed description would require a 
discussion of generators and relations for C*-algebras, a theory yet to be developed for real 



[Ui 



in M k (C(X)) with u' 1 



u 



u we find formulas for the 



C*-aglebras. Given u 
component function 

Uij = Uij{Xi, . . . , Xd+l) 

so that we interpret via functional calculus the meaning of 

(ax, . . . , ad+i) 

formulas for ^-representation a\, . . . , ad+\ of the ci-sphere in a C*-algebra. When 5 is small, 
we will find 

U = U{a x , . . . ,a d+ i) 

satisfies U^ 1 ~ U* = U and so have a well defined element of Kq(A). In particular, we were 
working in the GUE case over the sphere with 



u 



X 3 X\ - %x 2 
x\ + ix 2 —X3 



e M 2 (C(5 2 )). 



Moving to the 3-sphere, we have a Ax-element determined by the unitary 

G M 2 (C(S 3 )) 



u 



X3 + 1x4 
xi + ix 2 



—X\ + IX 2 

X3 — 1X4 
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and given a%, . . . , CI4 in a C*-algebra A that are a soft-representation of the 3-sphere we get 
a representative of a i^i-class in A using 

03 + —di + ia,2 

di + ICL2 Cf-3 — ICL4 

This leads to nothing of substance when A = M n (C), as this has K\ equal to zero. This 
pushing forward become interesting when we switch to the study of real C*-algebras. 

8.5. Real C*-Algebras. There are two competing mathematical objects forming the sub- 
ject called real C*-algebras. The one we fill emphasize are what we call C* T -algebras, which 
are C*-algebras with an extra operation. 

A C* T -algebra is a pair (A, r) where A is a C*-algebra and r is a *-isomorphism 

t : A — A op 

subject to the axiom r(r(a)) = a. We regard A op as having the same underlying set as A. 
The opposite of a C*-algebra is defined as having the original algebraic operations, except 
that the new multiplication applies the original multiplication in reverse order. We may 
just as well treat r as a unary operation a h-> a T that satisfies all the properties of being a 
*-homomorphism, except the multiplicative rule becomes 

(ab) T = b T a T , 

and with the axiom 

(a T ) T = a. 

Since * and r commute, we write a* T where many authors write a. The advantage of our 
notation is that one C*-algebra can have t\ 7^ r 2 that make it a C* T -algebra and lead to 
different a operators. In the case of A = M 2 tv(C) there are two possible r operations, the 
dual and the transpose, one of which leads to T being the expected entry-wise conjugate 
operation, the other does not. 

Mathematicians generally work with the real part of A, denoted §l T (A) and defined as 

Ur(A) = {a e A\a T = a*}. 

This is a *-algebra over the scalar field M.. We take here a definition: an R* -algebra is a 
normed *-algebra over M. that equals §l T (A) for some C* r -algebra A. Two essential examples 
in physics are 

9£ T (M 2JV (C)) = M 27V (10 

and 

^ (M 2JV (C)) = M N (U) 

where H is the algebra of quaternions. The isomorphism is not terribly elucidating, so we 

do not discuss it. We advocate, however, calling the condition = the quaternion 

condition. It is easy to see the matrices that satisfy the quaternion condition are all those 

of the form 

' A 'J? 
B A 

Remark 8.2. The standard term for an i?*-algebra is a real C*-algebra. A C7* T -algebra is 
generally called a "real" C*-algebra, where the quotation marks are part of the name [23J, 
and the additional operation is emphasized is a t— > a* T . We think our notation is more fitting 
in this context, as the a v-)- a T operation satisfies the axioms of the transpose A \-t A T . 
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Invertibles used 


Neutral element 


One-dimensional example 


Classical form 


K (A) 




(:-.) 


K (C) = Z 


Sig 


KM) 


no restrictions 


1 


Kx(C) = 





TABLE 7. Complex K-theory. 







Neutral element 


One- dimensional example 


Classical form 


K T (A) 


=fc T * 

ry ry ' ry ry 

JU iAs ^ JU JU 


(:-'.) 




l) = Z 


Signature or trace 


Kl(A) 


x* — x^~ 






) = Z a 


Determinent 


Kl{A) 


rf ry T ry^~ ry 

JU JU *! JU iAs 


(-'.:) 


K 2 (M 


) = z 2 


Pfafflan 



TABLE 8. Real i^-theory. 



An essential feature of complex C*-algebras is the matrices over them are again C*- 
algebras, in a unique way. In the case of C* T -algebras, we lose uniqueness, but have a 
canonical structure. Given {A,t) the r operation, also denoted r, of M n (A) is T <g> r on 
M n (R) O A = M n (A), so 

(a jk ) T = (a T kj ) . 

8.6. Kq K\ and K<i for real C*-algebras. For present purposes, we need to define K n 
for a C* T -algebra for n = 0, 1, 2 and n = 4,5, 6. We have no need for K% and K-j and omit 
them. We deal with n = 0, 1, 2 in this subsection in a manner that makes them most look 
like signature, determinant and Pfaflian. 

For n < or n > 7 we rely on the order-8 Bott periodicity, and so take K_ 5 to be defined 
to equal K 3 . In the next subsection we show how to deal with the cases n = 4, 5, 6 by an old 
trick involving the quaternions. 

Consider how we constructed i^-theory in the complex case, in Table [71 Within the sets 
of invertible elements, we made a choice on what symmetry to require — none or self-adjoint. 
We needed a neutral matrix to know how to embed a set of smaller matrices into a larger 
set. We define K^(A), with alternate notation K n (^R. T (A)) or K n (A,r) to fit in context, by 
repeating what we did in the complex case, but with the symmetries and neutral elements 
as indicated in Table The required subsets of invertible matrices are the following: 

GL4°](A,t) = {xe GL n (A) \x* =x T , x* = x} 

GLW(A, r) = {xe GL n (A) \ x* = x T } 
GLf 1 (A, t) = {x G GL n (A) \ x* = x T , x* = -x } 

For K we require 

x 

1 = [x] . 

-1 
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and only use elements of GL$(A,t) for n even. For K\ we require 

x i = N ■ 

These are standard descriptions of K and K\ of the i?*-algebra 

U T {A) = {a e A\a T = a*}. 

For K2 we require 



x 



1 
-1 



[*]. 



and only use elements of GL[^(A, r) for n even. This is not a standard picture of K2, so we 
must prove it is equivalent to the standard picture where the symmetries required on the 
invertible elements of M 2n (A) are 



T * 2 



and the same neutral element. This is as described in [26]. Notice that it is a feature of 
Bott periodicity that we avoid higher homotopy groups by changing symmetries. Even the 
first homotopy group introduces loops, which are bad computationally as they introduce an 
infinite object even when studying finite-dimensional matrices. 



Lemma 8.3. If u in a unitary in a C* -algebra then 

.2 



u 



-u 



u 



-1. 



Lemma 8.4. Suppose x is an invertible element in a C* T -algebra and 



JU Ju « Jb ,Xj • 



If x has polar decomposition x = up then xt = up' satisfies 



T * * 



for all t between and 1. 

Proof. With s = we have 

{up')* = (x(x*x) s )* = (x*x) s x* = -(xx*) s x = -x(x*x) s = - (up*)* . 

Working with polynomial approximations, one shows that for / applied to a normal element 
y we have 

(f(y)T=7(y*) 

and 

(f(y)Y = f(y T )- 

Applied here, the second tells us 

(up') T = {x{x*x) s ) T = {{x*x) s ) T x T = {{x T x* T ) s )x T = ((x*x) s )x* = (up')* 



□ 
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Lemma 8.5. Suppose x is an invertible element in a C* T -algebra and 

X = X , X = —1. 
// x has polar decomposition x = up then x t = up 1 satisfies 

T * 2 1 

x t ~ x t i x t ~ 1 

for all t between and 1. 

Proof. Again up 1 = x(x*x) s . From 

(x*x) (xx*) = 1 

we derive 

* / * \ — i 

XX = [X X) 

Therefore 

up t up t = x(x*x) s x(x*x) s = (xx*) s xx(x*x) s = —(xx*) s (x*x) s = —1. 

The fact that x\ = x* t is as before. □ 

Theorem 8.6. Let Abe a unital t-C* -algebra. The group (A) is isomorphic to the usual 
K 2 of the underlying R* -algebra 

A = {a E A\a T = a*} . 
The natural isomorphism is induced by a homotopy equivalence of 

(8.3) {x e GL 2n (A) | x* = x T , x* = -x } 
with 

(8.4) {x E GL 2n {A) | x* = x T , x 2 = -1 } 
that sends, in either direction, the class of x to the class o/ polar(rc) . 

Proof. The definition of ^2(^0) is 7Ti of the union of the GL^(A,r), for all n. Part of Bott 
periodicity, as in [26], states that there is a natural isomorphism of i^^o) with 7Tq of the 
union over even m of the sets in eq. 18.41 

The proof of the stated homotopy equivalence is contained in the previous three lemmas. 
We leave to the reader to verify that [x] [polar(x)] respects the inclusions as n increases, 
and is natural. □ 

8.7. K4 K§ and K§ for real C*-algebras. The key isomorphism in this subsections is, at 
the i?*-algebra level, 

K n (A x H) = K n+i (A). 

We will work, however, in the form of r-operations. This should be more familiar to physi- 
cists. 

We define K£(A), for 4 < n < 8 by 

K T n m 4 (A ®M 2 (C)). 

In the alternate notation, which is perhaps more clear in this instance, we take as definition 

K n (A, t) = K n -±(A <g> M 2 (C), t ® ft). 
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8.8. Examples. For a compact Hausdorff space X the C*-algebra C(X) = C(X, C) is equal 
to its opposite. We could take a variety of order-2 homeomorphisms of X to create a variety 
of r operations on C(X), but the only one with obvious relevance to the physics in this 
paper is the trivial choice, so r = id. The r-C*-algebra (C(X),id) has associated i?*-algebra 
C(X, E). 

The interesting element in the X-theory of C(S , M) is lives in = K^-d- For the 
2-sphere, we consider a i^-theory class [b] in 

K 6 (C(S 2 , R)) = K 2 (C(S 2 ) ® M 2 (C), id <g> Jt). 

Let fi(x, y, z) = x, f 2 (x, y,z) = y and /^(x, y, z) = z, where we regard S 2 as the unit ball in 
M. 3 . Using the Pauli spin matrices we define 

b = yj f r <g> ia r . 

Since the a T anti-commute, square to —1, have a\ = —a r and a\ = a r we discover 

where r = id <8> tt- Notice b is unitary, and has the correct symmetries to define an element 

[b] e K 2 (C(S 2 ) ®M 2 (C), id (8>U) 
but we cannot use 6 as it is in GL^ of C(S 2 ) <8> M 2 (C). Instead we must use 

b 



MJcr 



l(Jr 



There are topological ways to see [^©icrj;] is not the trivial element in this group, but it also 
follows from our example in |T], since triviality of [^©zcr^.] would have forced our Pfaffian-Bott 
index to be zero for triples of self-dual Hermitian matrices with small commutators. 

The finite dimensional example, that was essential in our study of 2D TR invariant systems, 

is 

K 6 (M 2N (C), ti) = K 2 (M 2N (C) <g> M 2 (C), (J ® ft) . 
Using the isomorphism above, and some standard facts in real .fT-theory [27J, we find 

K 2 (M 4JV (C), (-) T ) = K 2 (M 4 (C), (-) T ) ^ K 2 (M 4 (R)) = Z/2. 

We need more than just an example of an i^ 2 -class here. We need an explicit isomorphism 

iT 2 (M 4W (C),(-) T ) ->{±1}. 

That is, given an invertible, skew-Hermitian and skew-symmetric matrix, we must determine 
which of the two classes in K 2 it is in. The mapping is 

[X]^sgn(Pf(X)) 



The isomorphism 

then is 

where 



K 2 (M 2JV (C) ®M 2 (C),Jt® ft) -»- {±1} 
[X] ^ sgn (Pi(U^XU)) 

U = -^=(I + Z®a 2 ), 



and where X is invertible with 

x m = x t = _ x . 
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These Pfaffians can be effectively computed using the following theorem. Notice skew- 
Hermitian plus skew-symmetric implies the matrix is real. 



Theorem 8.7. If A in M 2 at(IR) is invertible and skew- symmetric, then there is an orthogonal 
matrix U and positive real numbers a±, . . . , a at so that det(£7) = 1 and A = UDU^with 



D 





-eai 










a 2 
a 2 












ax 
a N 



where e = ±1. Moreover 

Pf (X) = ea x a 2 • • • ajv- 
The eigenvalues of X are ±ai, ±a 2 , . . . , ±ajy with associated eigenvectors 

T 



u 







75 7§ 







Proof. The Schur factorization of real matrices gives us an orthogonal matrix U and a block 
upper-triangular matrix D, with 1-by-l or 2-by-2 blocks on the diagonal, and A = UDU T . 
Possibly by altering a single row in U, and a single column and row in D we can assure 
det(U) = 1. Since U is orthogonal, D is real, invertible and skew-symmetric, so B is actually 
block diagonal and no 1-by-l blocks and has only 2-by-2 blocks of the form 

-dj 

a 3 

We are mostly finished. The signs of the aj are arbitrary, but conjugating by a diagonal 
matrix with ±1 on the diagonal, and with an even number of negative elements, we can get 
the signs as indicated. 

The rest of the statements are standard, using in particular the formulas 

Pi(YXY T ) = det(F)Pf(X) 



and 



Pf 



Pf (X) = 



z 
-z 



= z 





' x 1 





( 





x 2 



^ = Pf ( X 1 ) Pf ( x 2 ) 



which holds for all Y and all skew-symmetric X, X\ and X 2 . Taken together, these formulas 
are a good definition of the Pfafnan. □ 
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8.9. Pushing forward real i^-theory. We need a version of Theorem 18.11 to apply in the 
GSE and GOE cases, to show that when we start with commuting self-dual Hermitians, 
the Z 2 index is trivial. As our numerical studies are for 2D and 3D GSE lattices, we state 
as plainly as possible a result that covers those two cases. We state for future reference a 
theorem to cover the GSE and GOE in all dimensions in as brief a form as possible. 

The maps between C* r -algebras that induce homomorphism of .fT-theory groups are the 
*-r-homomorphisms, meaning the algebra homomorphisms that preserve both the * and r 
operations. 

Theorem 8.8. Suppose X is a pathwise connected compact subspace of M. d . For any unital 
* -homomorphism 

cp : C(X) -+ M 2N (C) 

such that 

¥>(/) = &(f)f 

the induced homorphisms 



and 



are the trivial maps. 



^ : K 5 {C{X),id) -> K 5 {M 2N {C)4) 
<p* : K 6 (C(X),id) -> K 6 (M 2N (C)4) 



Proof. Using the structured spectral theorem, applied to <p(h r ) for h r the various coordinate 
functions of Euclidean space, restricted to X, we find a symplectic unitary U and points 
X\, . . . ,Xn so that 



<f(f) = u 



f(x N ) 



IP. 



All symplectic unitary matrices are homotopic, and all points in X connected by paths, so 
we may assume U = I and Xj — X\ and so 

¥>(/) = /(xi)/. 

If we let ipu denote the inclusion C M 2 at(C) and let ip : C(X) —> C denote the map 
/ 1 — y f{x\) then we have shown that, as far as i^-theory is concerned, ip is the composition 
loi/;, i.e. = i* o fa. Since 



tf 6 (C,id) = K 5 



ATi(H) = 



and 



K e (C,id) = K 6 (M) = K 2 (M) =0 
(see [27], for example) we see that ip* is the trivial map. 



□ 
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Theorem 8.9. Suppose X is a pathwise connected compact subspace ofR d . For any unital 
*-homomorphism 

V : C(X) -> M n (C) 

such that 

<p(f) = WW 

the induced homomorphism 

ip* : K q {C{X)M) ->• K q (M n (C),r)) 
must send to the trivial element the kernel of the map 

5 x :K q (C(X),id)^K q (C,-) 
where 5 X sends f in C(X) to f(x), for some chosen base point x in X. 

Proof. The proof is essentially as above, except that if r is the transpose, then we use a 
simultaneous diagonalization of commuting symmetric matrices by a real orthogonal matrix. 

□ 

We now return to the discussion of pushing forward K-theory by "almost homomorphisms." 
Consider once more 



u = 



X3 + ixi —xi + ix2 
X\ + 1x2 X3 — ixi 



GM 2 (C(5 3 )) = C(S 3 )® M 2 (C) 



but now notice it has symmetry 
This is defining an element of 



u = u 



K^CiS 3 , id) ® (M 2 (C), (t)) = K_ 3 (C(S 3 , id)) 



Now when we form 

U = U(a 1 ,a 2 ,a 3 ,a A ) = 



a 3 + ia 4 —ai + ia 2 
di + %a 2 a 3 — ici4 



e A® M 2 (C) 



for a 5- representation of the 3-sphere in (A, r), we obtain U that is invertible and has 

JJT®$ _ JJ* 

We have our induced element in 

K 1 (A^(M 2 (C)4))^K^(A). 



9. Numerical Implementation 

The main computations needed for the 2D and 3D GSE studies involve matrix functions 
of self-dual matrices and either the Pfafnan or determinant of a real matrix. The Pfafnan 
and most of the matrix functions depend on a structured factorization so we consider that 
topic first. Next we consider Newton's method for computing the polar of a matrix, and 
finally matrix functions of Hermitian and unitary matrices. 
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9.1. Factorization of self-dual matrices . Taking advantage of the block structure a 
self-dual matrix can create a block-diagonal matrix and so reduce the complexity of further 
computations. The algorithm we use is known as the Paige/van Loan algorithm. 

Theorem 9.1. (Paige/van Loan) If X is a self-dual, 2N-by-2N matrix, there is an order 
N 3 algorithm that produces a symplectic unitary U so that 

A C 



U^XU = 



A T 



for some A and C in Mat(C) with A in upper Hessenberg form and real numbers on the 
sub-diagonal. 



X 



Proof. We know 

" A C 
B A T 

with B T = —B. A standard fact in numerical linear algebra is that there is a Householder 
unitary V so that VB is a matrix with zeros in the first column below the first subdiagonal 
and such that V fixes e^ It follows that B' = VBV T also has zeros in the first column below 
the subdiagonal, so we consider the symplectic conjugation 



" V 


" 




' A 


C 




' V T 







" A' 


a ' 





V 




B 


A T 







vt _ 




_ B' 


A' T _ 



This unitary is symplectic, so we retain the self-dual structure and so B' is skew-symmetric, 
hence zero on the diagonal and so the first column of B' will be a scalar multiply of e 2 . 
Again by standard techniques we find a and f3 in C so that a 2 + (3 2 = 1 and 

W = a |e 2 ) (e 2 | + /3 |e 2 ) (e N+2 \ - P lejv+2) (e 2 | + a \e N+2 ) (e N+2 \ 
will be a unitary (a so-called Givens rotation) so that 



W 



A' 
B' 



a 



A" 
B" 



C" 
A" T 



will produce B" with the entire first column equal to zero. The construction of W is such 
that it will be a symplectic unitary so B" will be skew-symmetric and have first row zero 
as well. Another Householder unitary Q can be found so that Q fixes ei and QA" has first 
column with a real scalar in the subdiagonal and zeros below that. This fact that Qei = ex 
implies that QB"Q^ will also have first row and column equal to zero, so 



" Q 


" 




' A" 


C" ' 











- A m 


C" ' 





Q . 




B" 


A" T _ 







Q T 




B'" 


A"' T 



will have B'" with zeros in the first row and column and A'" having zeros in the first column 
below the first subdiagonal and a real scalar on the first subdiagonal. We use the symplectic 
unitary 

to find XJ\XU X is the row 1 and column + 1 in the correct form and JJ\ fixes both e x and 
ejv+i- This last condition allows for an iterative solution. □ 

If we apply the Paige/van Loan algorithm to a unitary or Hermitian matrix, we get an 
even better outcome, a block-diagonal matrix. 



Q 


" 


w 


' V 


" 





Q . 





V 
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Theorem 9.2. If 



X 



is self-dual and normal, then A is normal, D = 
Proof. Since X is self dual, D = A T . Therefore 

A f 



A C 
D 

= A T and C 



0. 



and 



X*X 



xx ] 



C* A 

A C 
A T 



A C 
A T 

A f 0_ 
Ct A 



' At A _ A^C 

AA T + C*C 

AA^ + CC^ CA 
A T C*t A T A 



The normality of X tells us AA* < A^A. For matrices, hyponormal implies normal so At A = 
AA^ and &C = 0. Thus C = 0. □ 

We have now a method of diagonalizing a unitary or Hermitian self-dual matrix X. From 
the Paige/van Loan algorithm, we obtain a symplectic unitary U so that 



X = U 



Y 

o r T 



c/t 



with y being unitary or Hermitian. In the unitary case use use standard algorithms to obtain 
a Schur factorization Y = VDV^ with D upper triangular and V unitary. As D is unitary, 
it is also diagonal. In the Hermitian case, we use a standard eigensolver to diagonalize Y. In 
either case, 

" D 
D 



X = U 



V 0_ 
V 

so the desired symplectic unitary is 



rt 
o 



o 



c/t 



u 



V 0_ 
V 



In practice, when we diagonalize a self-dual self-adjoint matrix, specifically a Hamilton- 
ian in a GSE system, we store only the left half of the all intermediate matrices that are 
either Hermitian self-dual or symplectic unitary. This allows the structured eigensolver to 
be competitive in speed with unstructured eigensolvers. The timing of such algorithms is 
architecture dependent, and very sensitive to the test matrices, but Table [9] gives some ideal 
of accuracy and speed. 

We need to deal with approximate unitary matrices as well, so have the following variation. 



Theorem 9.3. // 



X 



A C 
D 



is self-dual and | JTtX 
Proof. Since X is self dual, D 

X*X-I 



< 5 then D = A T and \\A ] A 
A T , and 



/ < 5 and ||C|| < V2S 



"At " 




"A C 




" AtA - 1 


A^C 


Ct A 




A T 




0A 


AA T + CtC - / 
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AT 
IV 


Time (sec.) 


^ H-UDW 


max — Xj\ 


1 IT TTDTTi 
YgT5- H — UL>U 1 




stand. 


struct. 


stand. 


struct. 


stand. 


struct. 


adjusted standard 


500 


7.6200 


6.4800 


0.0176 


0.0129 


0.0066 


0.0000 


0.7516 


1000 


62.530 


60.020 


0.0259 


0.0206 


0.0112 


0.0000 


2.0157 


1500 


207.68 


200.17 


0.0307 


0.0225 


0.0115 


0.0000 


6.0742 


2000 


488.13 


472.22 


0.0360 


0.0299 


0.0162 


0.0000 


10.796 



Table 9. Comparison of eigensolvers for W = H = in M2at(C). In the 
standard case, using ZHEEV, D is diagonal, and U is unitary. In the structured 
case, U is a symplectic unitary and D is diagonal with doubled eigenvalues. 
The last column reflects the error arising from adjusting eigenvector pairs in 
the output of ZHEED to force the expected symmetry. Each reported average 
is over 10 test matrices. 



which implies ||^^4 — /|| < S. As we are working in finite dimensions, we know \\AA^ — l\\ < 
5 and WXX^-lW < 5. From 



XX ] -I 



' A 


C 




' At 


" 







A T 






A 


- I = 



we derive 



so 



\C\ 



AA^ + CCt _ / CA 
A T & A T A - I 



\\AA* + CC* -I\\ < 5 
|CC t || < WAA 1 + CC ] - III + \\AA ] - ill < 25. 



□ 



9.2. Factorization of skew-symmetric matrices. We need to calculate the Pfaffian of a 
Hermitian, skew-symmetric matrix Y in M 2 tv(C). Rather than factor Y as in Theorem 18.71 
we set X = —iY, we compute the Pfaffian of the real matrix X and use formula 

Pf (iX) = i N Pf (X) X G M 2N {C). 

For a real, skew-symetric matrix, we can combine standard algorithms to get a factorization 
that exposes the spectrum and Pfaffian at the same time. 

Theorem 9.4. There is an order N 3 algorithm which, for X in M 2 at(1R) that is skew- 
symmetric, calculates Pf(X), an orthogonal matrix U and real numbers ai, . . . so that 
X = UDU T with 



D 



a\ 
-a\ 






a 2 
-a 2 













-a N 



a N 
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Proof. We start with a Hessenberg decomposition X = QYQ T with Q orthogonal and Y 
having all zeros below the first subdiagonal. As Y = Q T XQ must be real with Y T = —Y, 
we in fact have Y is tridiagonal, indeed of the precise form 



(9.1) 



Y 







Cl 



-6i 



6i 

c 2 
-c 2 b 2 

b 2 •■• 





-b N -i 

-Cat 



Cat 





The Hessenberg decomposition can computed in LAPACK by DGEHRD. 

We now follow ideas of Golub [37]. Let P be the permutation matrix corresponding to the 
shuffle so that 



P T YP 



C 
-C T 



where 



C 



Cl 

-h c 2 



CjV-l 



Take a singular value decomposition C = WEV with W and V orthogonal and 

«i 



«2 



Q>N—1 



a N _ 



for Oj > 0. Then 










c ' 




" w 




E ' 


V 















-E 



and so 



P 



V 



c 
-C T 



p T = p 



E 



-E 
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where 



D 



ai 
-at 

a 2 

-02 

'•. 




a N 
—a N 



Next 



so we are done with 



P 



V 



P T Q T XQP 



W 



V s 



P T = D 



U = QP 



W 



V' 1 



P T . 



□ 



9.3. Stuctured polar decomposition. As argued by Higham [36], a fast and accurate 
way to compute the polar of a matrix is via Newton's method. Suppose X is an invertible 
complex matrix. To compute U unitary and P positive with X = UP we set X\ — X and 



X k +i - - (x n + [x n x ) f 



and get U = lim^Xfo. The convergence is quadratic in the operator norm, in in practice we 
used three iterations. 

Theorem 9.5. If X 9 = ±X then X[ = ±X k for all k and 

polar (Xf = ±X. 
IfX* = ±Xt then X[ = ±X ] k for all k and 

polar (Xf = ±X ] . 
IfX T = ±X then X^ = ±X k for all k and 

polar (X) T = ±X. 
IfX T = ±Xt then X T = ±X ] k for all k and 

polar (X) T = ±X f . 

Proof. From the identities (X^ 1 = (X" 1 )* and (X») t = (X^f we deduce 



X 



fc+i 



- xl 



X 



-1 



and 

Similarly 



polar(X) tt = polar(X a ). 
^i = i(^f+(W)" 1 ) t ) 
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and 



From (Xt) 1 



polar(X) T = polar(X T ). 
(X- 1 ) t and (X^ = X we deduce 



X l+i 




and 
Clearly 



polar(X) t = polar (X f ). 
-X k+1 = - (-X n - (X" 1 ) 1 ^ 



and so 

polar(— X) = — polar (X). 
All eight versions of the theorem follow from some combination of these formulas. 



□ 



The most accurate way to find the polar of a self-dual matrix X is to apply Newton's 
method to X. A faster, less accurate method, which we utilized, is to apply Paige/van Loan 
to get symplectic unitary Q with 

A C 
A T 



X = Q 



Q 



apply Newton's method to A to find U unitary with U m polar(A) and then use 



Q 



U 
U T 



which will be a self-dual unitary, close to the polar of A when X*X I. 

We can take logarithm of U by using a Schur factorization, as we discuss in the next 
subsections. This produces for us a Hermitian matrix Y with e 2mY ~ U, and from here form 



G = Q 



Y 
F T 



Q 



which is Hermitian self-dual with e 2mG ~ X. We summarize our ability to find self-dual log- 
arithms of self-dual approximate unitary matrices in Table, using test matrices that are very 
close to being unitary. We compare our results to what results from using an unstructured 
Schur factorization. Clearly we could do better in calculating the log in the standard case, 
by attempting to pair the eigenvalues and adjust the eigenvectors before taking log. The 
alarming errors are a warning that something must be done when applying discontinuous 
functions to normal self-dual matrices to get any semblance of a self-dual matrix out. 

9.4. Matrix functions of structured Hermitian matrices. Given a measurable function 
/ : C — > K. and a normal matrix X the definition of the matrix /(X) tells us to compute it 
from any unitary diagonalization of X, 



Ai 



/(Ai) 



(9.2) 



X = U 



f(X) = u 



c/t. 



A r 



/(A n ) 
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A T 

Jy 


i T 7"t T T T 

WU-I 


Time (sec.) 


U-QDQ^ 




rrti TT 




samples are exactly 
self- dual 


stand. 


struct. 


stand. 


struct. 




etrn of 

Oil \JLKj u • 


500 


0.1028 


25.002 


16.799 


0.0760 


0.0743 


2.8309 


3.3329 x 10" 14 


1000 


0.2457 


183.11 


140.52 


0.1674 


0.1754 


2.8049 


5.1817 x 10" 14 


1500 


0.4298 


641.89 


465.04 


0.2876 


0.3079 


2.8803 


5.5677 x 10" 14 


2000 


0.6618 


1480.1 


1086.8 


0.4381 


0.4743 


2.8320 


6.5637 x 10~ 14 



Table 10. Comparison of eigensolvers/logarithm for W = U~ x , U* = U^ 1 in 
M 2 at(C). In the standard case, using a Schur decomposition of U found with 
ZGEES, D is a unitary diagonal, and Q is unitary. In the structured case, 
Q is a symplectic unitary and D is self-dual unitary diagonal. In each case 
H — Qt-L i g Q Each reported average is over 10 test matrices. 



Notice f{X) is Hermitian since we asked that / be real- valued. We defined P this way 
from the Hamiltonian, using the function X(-oq,e f ], and saw in that case that P will be self- 
dual when the Hamiltonian is self-dual. This sort of preservation of structure is true more 
generally. When / is not continuous we will generally need to numerically compute f(X) 
using a structured unititary factorization, as in the proof of the following. 

Theorem 9.6. Suppose f : C — > R is measurable function X is a normal matrix X. If 
X e M 2N (C) is self-dual then f(X) is self-dual. If X T = X then f(X) T = f(X) and so 
f{X) is real symmetric. 

Proof. When X is self-dual, we can find a symplectic unitary U so that the second half of 
the eigenvalues repeats the first. This implies that 

/(Ai) 



f{X) = U 



/(A 



N 



/(Ai 



/(A 



U* 



N) 

is self-dual. An algorithm to find U is discussed in §9.11 

Notice that when X T = X and X is Hermitian, we are assuming X is a real matrix. In 
this special case, we can use standard algorithems to find U in f)9.2p that is real orthogonal. 
The fact that f(X) is then real symmetric is now obvious. The more general form of the 
second claim will arise when applying the function /, g and h as in Section .... to symmetric 
unitary matrices, as wil be necessary in the study of GOE systems on the torus in dimensions 
4, 6 and 7. This is a theoretical issue at this time, so we don't worry about an algorithm. 
We write X = A + iB for A = \{X* + X) and B = f(Xt - X) and find that A and B 
commute (since X is normal) and are real symmetric since X T = X. There is then a real 
orthogonal matrix U that diagonalizes both, and so diagonalizes X. We again can use (19.21) 
to show f(X) is real symmetric. □ 

9.5. Matrix functions of almost unitary matrices. The definition of f{X) corresponds 
to more naive notions of a function of a matrix when / is a polynomial, or even a Laurent 



(:>(> 
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function, i.e. a polynomial in X and X 1 . Of course the latter will require X be invertible. 
For example, if 

f{x) = 2x 2 + x- 1 

then 

since given (19. 2p we have 



f(X) = 2X 2 

K 1 



x- 1 



x- 1 = u 



A; 



and 



r \2 



x 2 = u 



A? 



\2 



Z7+. 



Unless high powers are involved, for such functions the diagonalization will be slower than 
those more naive calculations. More importantly, when X is approximately unitary, we have 

x^x- 1 . 

Define X™ = X n for n nonnegative, and X( n ^ = (Jft)". As long as \k\ is small, we have 
reasonable approximations X^ « X n and 

rn 

fix) « ^ a w 

fc=— m 

when m is small and / is the Laurent polynomial 

m 

m = J2 



CLyi X 



The functions we used when mapping the torus to the sphere were 

/( A ) = lis sin(arg(2;)) + -f, sin(3 arg(z)) + ^ sin(5 arg(z)) 



9W 



h(z) 



128 

arg(z) G [t, 



1 - {fi z )f aT S( z ) 



I 31 
4' 4J 

[I 31 

U' 4J 







arg(2;) ^ [ x :i 



4' 4J 



where now we are using a complex argument A or modulus one, so that A (n) = A n . These 
functions were selected to be close to degree-five Laurent polynomials. Indeed / is exactly 
and order-five Laurent polynomial, while estimated the Fourier series of g and h to arrive at 

555 



m = £ 



a n z {n \ g(z) 



Kz {n \ h(z) « J2 c " z 



(n) 



n=— 5 n=— 5 n=— 5 

with the coefficients as in Table [TT1 
Rather than computing 

/(polar(£/)), #(polar([/)), %olar([/)) 
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n 





1 


2 


3 


4 


5 


CL n 





-150i 
256 





256 





-3i 
256 


b n 


0.20205 


-0.17994 


0.125655 


-0.06601 


0.023445 


-0.0038855 


Cn 


0.20205 


0.17994 


0.125655 


0.06601 


0.023445 


0.0038855 



Table 11. The coefficients denning /, g, and h, extended to negative indices 
by a_ n = and b_ n = b n and c_ n = c n . 



we compute 

5 5 5 

F = a nU {n \ G=J2 h nU [n \ H=J2 C nU [n) 

n=— 5 n=— 5 n=— 5 

to save time. We give up exactness in expected relations, so get F 2 + G 2 + H 2 / and 
G# « 0, but when \\WU - l\\ is small, the errors \\F — /(polar(i7))|| are small so we will 
not change any index. Computing a fifth power of a matrix introduces little error, so when 
U is self-dual within machine precision, so will be F, G and H. 

10. Discussion 

We have given a general procedure for obtaining topological invariants of free fermion 
systems by mapping to a problem in C7*-algebra. The general approach for systems in the 3 
classical universality classes on a sphere is to form a projector P onto the space of occupied 
states, and then to project the coordinate matrices of the system into the occupied band. 
These matrices approximately commute and the sum of their squares is approximately equal 
to the identity. We regard these matrices as describing a soft sphere S d . We then look for 
topological obstructions to approximating these projected matrices by exactly commuting 
matrices. To compute these topological obstructions, we form the operator B in (12.11) . This 
operator B approximately squares to the identity. Thus, we can view this operator B as 
describing a soft zero dimensional sphere. In this paper we have described an approach to 
computing topological properties of B in any dimension and six of the ten symmetry classes. 

Previously, we implemented a special case of this procedure in two dimensions [1]. The 
algorithm based on C*-algebra allowed us to study much larger systems than techniques 
based on considering a flux torus [11] in two dimensions. 

In this paper, we have implemented this procedure numerically for three dimensional 
time-reversal invariant insulators, considering systems up to 12 3 sites, with 4 states per site, 
for a total of a 6912-dimensional Hilbert space. We have found numerical evidence for an 
unexpected transition in the thermodynamic limit. In the presence of disorder, the system 
has localized states close to the band edge and close to zero energy, with delocalized states 
at intermediate energies. The invariant, averaged over disorder, also appears to display 
transitions in the thermodynamic limit, being equal to +1 for all samples for Ep near the 
band edge and equal to —1 in all samples for Ep near zero, and it fluctuates from sample to 
sample over an intermediate range of energies. That is, there appear to be critical values of 
the energy, E c , E' c with E c < E' c < such that the invariant fluctuates from sample to sample 
for E c < E F < E' c but the invariant is not fluctuating for Ep < E c or E' c < Ef < 0. However, 
while the energy at which the localization transition occurs near zero energy appears to 
coincide with the energy E' c of the transition in the invariant, the transition at energy E' c 
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appears to be unrelated to the localization transition. If so, this reveals an expected phase 
transition which has no signature in the localization properties. 

Topological properties of the soft zero dimensional sphere are described by one of three 
properties of the operator B. These are the Pfaffian, the determinant, and the number of 
positive eigenvalues of B. In one specific example, in the two dimensional time- reversal 
invariant case, we applied a similarly transformation to -B[TJ[9] to make B anti-symmetric, 
and then computed the Pfaffian. In the three dimensional case, we chose a particular set of 
7 matrices to make B block off-diagonal, and then computed the determinant of one block 
of B. 

For systems in the chiral universality classes, Wannier functions are no longer the defining 
feature of whether or not a given Hamiltonian can be connected to a trivial Hamiltonian. 
Instead, we have a unitary (representing the upper-right of the spectrally flattened Hamilton- 
ian) which approximately commutes with the position matrices. This alows us to construct 
invariants for 6 out of the 10 universality classes from invariants in C*-algebra. Some of 
these invariants were known invariants in C*-algebra, such as the invariants in the GUE 
case, while other invariants, such as the invariants of self-dual matrices in two [9] and three 
dimensions were new. We leave the question of the remaining 4 universality classes for the 
future, but we expect that these problems also have interpretations in terms of C*-algebra 
which will also lead to new algorithms.. 
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